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SECTION I 

INTRODUCTION AND BAOCGROOND 


This user's guide, describing the Non- Intrusive Parameter Analysis 
Software (NAS), is intended as an introductory text. Its organization 
comprises an outline of topics for software operation by a new user. Some 
features are only partially presented at fiist; their complete explanation 
is postponed to the point where additional Information will not be con- 
fusing or out of context • This guide Is oriented toward a user who is 
familiar with FORTRAN. Familiarity with FORTRAN will be required to exer- 
cise some of the special features of this software; however, a knowledge 
of FORTRAN is not required for basic program operation. It is assumed 

that the user is familiar with the operating procedures for the computer 
system in which this software resides. 

NAS incorporates many features over which the user has direct control 
and allows the user to adapt the software package to the identification 
problem at hand. The capabilities of this software include the ability to 
identify linear as well as nonlinear relations between input and output 

parameters; the only restriction is that the input/output relation be 
linear with respect to the unknown coefficients of the estimation equa- 
tion. The output of the identification algorithm can be specified to be 
in either ^ time domain (I.e., the estimation equation coefficients) or 
in the frequency domain (i.e., a frequency response of the estimation 

equation). The frame length ("window") over which the identification 
procedure is to take place can be specified to be any portion of the input 
time history thereby allowing the freedom to start and stop the identifi- 
cation procedure within a time history. There also is an option which 
allows a "sliding window." Under this option the window is "slid" to 
obtain new data, while the data at the beginning of the window is 

"dropped" as the new data is added. This option gives essentially a mov- 
ing average over the time history. 
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The software package Includes a user-defined subroutine which can 
manipulate the Input data of the Identification algorithm Into a conven- 
ient form. Another user-defined subroutine specifies the form of the 
assumed solution. Using this subroutine It is possible to run many as- 
sumed solutions simultaneously, which is advantageous when the exact form 
of the solution is unknown. 

ORGANIZATION 


Section II describes the fundamental concepts of the Non- Intrusive 
Parameter Identification Procedure (NIPIP). It outlines what Is required 
of the user to understand and implement this technique. Section I T also 
provides a discussion of the special options of the software that add 
scope and convenience to the user's ability to manipulate and Interact 
with the program. As the user becomes more familiar with the concepts of 
the NAS software package, this section provides several additional fea- 
tures that are of immediate usefulness. 

Sections III and IV comprise the majority of the information for the 
operation of the software package. Section III defines the program var- 
iables with which the user must be familiar to understand and operate the 
software. This section contains a description of the program's control 
flags and their use in routing the program during execution. It also 
contains a description of the input and output specifications and the 
internal warning messages which alert the user to possible error condi- 
tions and solutions of questionable validity. Section IV provides some 
illustrative examples in which problems are described from setting up the 
estimation equation through describing the output of the algorithm. Al- 
though this section does not examine every conceivable possibility, the 
user should be able to understand the operating procedure and run the 
software after completing these sections. 

Appendix A describes a small perturbation simulation of the F-8 air- 
craft which was used to generate some of the examples used In 
Section IV. Appendix B contains detailed Instructions for running the NAS 
program. Appendix C is a summary of the NAS program variables. Program 
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SECTION II 

BACKGROUND AND THEORY OP OPERATION 

A. DESCRIPTION OF NIPIP 

This section describes the basic concepts of the Non-Intrusive Param- 
eter Identification Procedure. The general procedure presented here was 
demonstrated in Ref. 1 and was extended to real-time computation in 
Refs. 2 and 3. This method involves the assumption of a general parameter 
identification strategy (a mathematical structural model with undetermined 
coefficients) and determines the coefficients of the strategy by matching 
the model to the generated data using a running least-squares estimation 
technique (i.e., a multiple linear regression). 

First , the general concepts of the least-squres approach are dis- 
cussed. Next the method by which a general form of an assumed pilot 
control strategy is stated. Finally, a description of the NIPIP algorithm 
is presented. All of this is in preparation for the presentation of the 
program variables and Input/output specifications which follow in the next 
section. 


1. Basic Cbncept 8 

The numerical analysis concept central to the procedure presented here 
consists of applying a running least-squares estimation technique. A 
specific example of least-squares estimation occurs in curve-fitting prob- 
lems where it is desired to obtain a functional form of some chosen order 
that best fits a given set of discrete measurements, Ref. 4. Thus the 
least-squares problem consists of correlating a dependent variable, y, of 
a given system with one or more elements of an independent variable vec- 
tor, That is, it is assumed that y and are related by an expression 
of the form 
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y * fjCj + f 2 c 2 + ... + fjCj ... + f a c m 



( 1 ) 


or y i £c 


( 2 ) 


where fj, fj, ••• are selected variables from the state vector jC or ex- 
plicit functions thereof, and the c j' s are constant coefficients relating 
y and f y In Eq* 2 JF is a row-vector of the f j, and £ is a column-vector 
of the Cj« y and can be composed of a control or state, past or present 
discrete values, and linear or nonlinear functions thereof, e*g*, 
”y * F c” could include the following forms: 


or 


or 


or 


e 

6 


e 1/2 


6c i 

§ c l + c 2 

^ C 1 + 0C 2 + ® c 3 + ^ c 4 + c 5 

2 n 5/2 

^1 + n c 2 + 7+JT c 3 + c 4 


If there exist sets of discrete measurements for y and _F, and the 
difference between the actual system output and the output of the assumed 
model given by Eq* 1 is given by v, then 


F c + v 
-n — n 
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Where c is the estimated value of c given by Eq* 1, and the subscript n is 

A 

used to denote each set of y and It Is assumed that there exists a £ 

which will relate all values of y n to the Fp. 

A 

”0 arrive at an estimate for £, a minimum of m sets of measurements 
are required (where m is the number of degrees of freedom in Eq. 1). The 

A 

solution for £ is found by adjoining these sets of discrete measurements 
is follows, 


- 


- 


“ * 

y l 

y 2 


Il 

*2 


v i 

V 2 

m 

m 

• 

C + 

• 

• 


• 

— 

• 

• 


• 


m 

1 

as 

1 


_In 


V N 


(3) 


where N is the number of measurements and N > m. Equation 3 can be writ- 
ten more concisely in a matrix notation as 


1 * H £ + v 


(4) 


where £, H, and v_ are defined by inspection. Note that H is an N by m 
matrix. 

A. 

The least-squares estimate for the unknown coefficients, £, is ob- 
tained by minim" ing the sum of the squares of £ where 


v * 


v - He 


(5) 
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The sum of squares of a vector is given by the vector inner product. Thus 
we wish to minimize the scalar J, where 


J 


" H £ ) T (£ - H £ ) 


( 6 ) 


Minimization of a scalar with respect to a vector is obtained when 


3J 

Jc 


0 


(7) 


and the determinant of the Hessian matrix of the second partial deriva- 
tives of J is positive, semidefinite 


> 0 ( 8 ) 



Differentiating J and setting the result equal to zero yields: 


H T Hc_ » H T 2_ 


(9) 


It is readily shown that the determinant of the Hessian matrix of J, with 
* 

respect to c_, is positive semidefinite, and thus Eq. 9 does, Indeed, de- 
fine a minimum. When H^h possesses an inverse (i.e., when it has a non- 
zero determinant) the least-squares estimate is 


£ - (H T H) _1 H T £ 


( 10 ) 
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Implied in the preceding development is that all available measure- 
ments are utilized together at one time (i.e., in a so-called ’’batch 
processing” scheme). 

Although the estimate of the unknown coefficients can be obtained 

using Fq. 10, it usually requires vast amounts of storage to save all of 

the accumulated data which is needed to form the H matrix. This storage 

requirement is sometimes prohibitive to the application of this scheme. 

However, it is easily shown that a summing of the measurements as they 

T T 

become available will yield the same H H and H matrices as those ob- 
tained by adjoining each set of measurements. Thus, 


and 


T 

H H 


N 

Z 

n=l 


T 

F F 
— n — n 



N 

Z 

n=l 



( 11 ) 


( 12 ) 


This gives the recursive least-squares estimate as 



N 

Z 

n=l 


T 

f y 

— n J n 


(13) 


where N is the number of data points up to the given point in time where 
the sum Is truncated. 

The following features of this solution can be noted. First, the main 

computational task consists of updating summations of products of the m by 

m matrix _F_F . The only storage requirement is that these summations alone 

be saved, not all of the accumulated data as with the batch processing 

scheme. Second, the most complex computational task is the inversion of 
^ T 

the Z matrix; however, this is still not a serious complication. 

n-1 
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2 m Assumption of a General Pilot 

Control Strategy Fora 

One use of this identification scheme is to identify the pilot's con- 
trol strategy by matching an assumed strategy to pilot-generated data, 
Basic and applied research going back more than two decades (Refs. 5 
through 10) has demonstrated that the human operator's control strategy 
can be fairly accurately described by linear differential equations. When 
the control task is time-varying, such as controlling an aircraft on final 
approach, then the coefficients of the differential equations will, in 
general, also be time varying. The major interests here, however, are the 
states that must be perceived and acted upon in order to perform the con- 
trol task. 

A general control task is depicted in Fig. 1. 

External 



Figure 1. An Example of a Pilot-Vehicle System 

This example represents the use of the controls, and the aircraft 
states, JC, for the identification of the pilot's actions. There are sev- 
eral concerns inherent in this kind of approach. The main concern is that 
since the raw data used in the recursive least-squares (i.e*, curve fit* 
ting) procedure consists of the controls, j5, and the aircraft states, _X, 
as shown in Fig. 1, there can be confusion as to whether the correlation 
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between j5 and _X reflects more accurately functions describing the pilot, 
Yp, or those describing the aircraft controlled element, Y Q . However, 
since Eq. 13 guarantees a minimum of the scalar, J (Eq. 6), and this mini- 
mum is unique such that the functions describing the controlling element 
and those describing the inverse of the controlled element do not yield 
the same minimum for J, then Eq, 13 will yield the solution for which J is 
minimized. If the assumed model form more appropriately describes the 
controlling element, the guarantee of minimizing J assures that the corre- 
lation obtained using Eq. 13 will identify the controlling element. 
However, this guarantee of minimizing J does not guarantee that the solu- 
tion will reflect either element if the assumed form of Eq. 1 does not 
adequately describe either element which is to be identified. It has been 
found that so long as the functional form of the assumed control law, 
i.e., y * c_ 9 is appropriate for the piloting technique being used, this 
procedure does identify the control function of the p*lot, Yp. This fact 
was demonstrated in Refs* 2 and 3 and will be shown ^.o be the case In a 
later section. 

The advantage of this procedure is that it is mathematically simple, 
as stated earlier, but its success depends .pon the judicious choice of an 
assumed pilot control strategy, Yp. This, in turn, depends upon the ex- 
perience of the experimenter, knowledge of the control task, and a 
thorough understanding of the dynamics of the controlled element, Y c . It 
is essential that the analyst choose a likely candidate control strategy 
in order for the subsequently estimated parameters of the control strategy 
to be valid. 

Some common examples of pilot control loop structures are shown in 
Figs. 2 and 3. These examples represent variations in the use of controls 
and variations in available displays or pilot-perceived aircraft states. 
The loop structures in Fig. 2 are examples of single input control-loop 
structures. Figure 3 gives examples of parallel loop structures 
(Fig. 3a), series loop structures (Fig. 3b), and control crossfeeds 
(Fig. 3c). The "Y" symbols within blocks in the loop structures represent 
describing functions (or control laws) which characterize the pilot's 
control strategy. It is important to recognize that each of these control 
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Stick 
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n. Altitude Rate of Change Task (Discrete Task) 


Vertical 



c. Landing Maneuver (Discrete Task) 



Figure . Three Common Control Loop Structures 
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loop structures can be handled and identified by the procedure being 
considered. 

The first step in applying the identification procedure is to formu- 
late the pilot control strategy for the defined control task. There is an 
element of engineering judgment in the choice of the structure of a mathe- 
matical function which appropriately describes the general input/output 
relation (I.e., the form of Yp) for the defined control task. Indeed, it 
is the analyst's exercise of a priori knowledge and experience which al- 
lows the judicious choice of the pilot control strategy. For example, a 
pure gain describing function is the most desirable so long as adequate 
pilot-vehicle performance can be obtained (Ref. 6)* Various combinations 
of lead or lag or other equalization can also be adopted by the pilot 
depending upon the response of the vehicle and the overall level of per- 
formance desired. 

The input/output characteristics which the pilot is likely to exhibit 
are readily expressed using Laplace transformations (transfer functions) 
of the differential equations describing the pilot control strategy; how- 
ever, this transformation yields a continuous model of the pilot. Since 
the identification procedure involves a digital computer, it is necessary 
to formulate the relationship, Yp, as a discrete relationship (i.e., a 
finite difference equation with undetermined coefficients). One way this 
is accomplished is by applying the z-transformation to the continuous 
system transfer function to obtain an equivalent discrete transfer func- 
tion. Table 1 presents some common z-transforms of continuous systems, 

along with the number of undetermined coefficients. This table is by no 
means complete and additional z-transforms can be found in Refs. 11 
and 12. 

3. NIPIP Algorithm 

Having discussed, in generalities, how to specify the form of the 
pilot control strategy, Yp, one can form the difference equation (estima- 
tion equation) and apply the running least-squares error estimation 
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TABLE 1 



Entry 

No. z-Transforra 


1 


2 


3 


4 


5 


6 


7 


8 


9 


1 - z 


1 - ajZ 


-1 -2 
b j z + b 2 z z 

” -2 
1 - a jZ - a 2 z 


bjZ * + b^z ^ 

_ — 

1 - a } z - a^z 


(z' P ) 


v -1 . v -2 

b.z + b„z 

V ^ 

1 - a^z 


v -1 . v -2 

bjZ + b 2 z 

-I -2 

1 - - a 2 z 


-1 -? -3 

b^ + b 2 z L + b ? z 

" Ti ^2 -3 

1 - ajZ - a 2 z - a^z 


. "I v ~2 , -3 , -4 

b.z + b„z + b_z +b,z 
1 2 3 4 

, -1 -2 -3 

1 - ajZ - a 2 z - a^z 
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algorithm in order to solve for the unknown coefficients of the estimation 
equation. In the following discussion of the NIPIP algorithm, a control 
task with an assumed pilot control strategy is used as an example. 

The control task is taken to be pitch attitude regulation (Fig. 2), 
with a lead/lag pilot control strategy (Entry 6 in Table 1). Thus, 


V 8) 


C + 1) 

K p ( TjS + 1 


exp(-i s) » 
e 


Pilot Series 
Gain Equilization 


Pure Time 
Delay 


6(b) 

9 (s) 
e 


(14) 


For this example, there is assumed to be no time delay, i.e., T g * 0. 

Once the control task and assumed control strategy have been speci- 
fied, the general form of the estimation equation can be determined. 
Using the 2 -transform given in Table 1, the discrete relationship for the 
lead/lag strategy is given by 


V z) 


. -1 , -2 
bjZ + z 


1 - a^z 


-1 


i(z) 

9 (z) 
e 


(15) 


The coefficients b Q , and b^ involve combinations of the parameters in 
the pilot's strategy, and it is not, in general, easy to use them directly 
to interpret the attributes of the pilot control strategy described in 
Eq. 14. For this reason it is often necessary to use the frequency domain 
to interpret the coefficients of Y p (z). That is, 


Y p (s«jo>) - Yj/z-e^ 1 ) 


(16) 


where w is frequency and T Is the discrete sample period. Effective val- 
ues of the Laplace transform frequency domain parameters K^, T^, Tj, and 
i e can, in turn, be approximated from the frequency response. 
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The general procedure for identifying the coefficients of the estima- 
tion equation is started by expressing the z-plane transfer function for 
Yp as the finite difference equation: 


or 


z _1 + b. 


e z _1 + 


V 


-2 


(17) 


a l 6 n-i + b l 0 e n-1 + b 2 0 e 


n-2 


(18) 


where the subscripts n, n-1, n-2, etc., denote the values of <5 and 9^ at 
times t, t-T, t-2T, etc., respectively, and T is the sampling period. 
Expressing Eq. 18 in the form of Eq. 2 gives 


y n 


(19) 


where y n = and the vector is described by 


In 


[«n-l 


■n-1 e n-2 


( 20 ) 


and the c vector is given by 


c 



( 21 ) 
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The values of the estimation equation coefficient, are determined using 
Eq. 13, 


N T , N 

c * ( Z F F ) I F 6 

— v , — n —n , — n n 

n*I n«l 


( 22 ) 


Equation 22 can be rewritten as 


where 


and 


c - (H T H] _1 H T <5 


(23) 


T 

H H = 


N T 
IFF 

i “i* 

n*=I 


N 

I 

n**I 


n-1 


n-1 

"n-2. 


[6 6 

1 n-1 i 


e ] 

®n-l e n-2 


(24) 


H T 6 
— n 


N 

l 

n«l 


F 6 
— n n 


N 

l 

n=l 


n-1 


n-1 

L e n-2j 


f 6 ) 

' n' 


(25) 


Tn future references, the notation of Eq. 23 will be used. The solution 
for c was shown earlier to minimize the squared error between the estl- 

A 

mated control, 3 = _F c_, and the actual control, 6. 

A measure of the goodness of the fit achieved by Eq. 13 can be com- 
puted as the integral square error between the actual and fitted data 
normalized with respect to the integral square of the actual data, 
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Z [« - « ] 2 
1 11 11 

R 2 - 1 - (26) 


where 6 * Hc_ is Che outpuC of Che model, N Is Che number of data points up 
Co Che given poinC in time, end R 2 is referred to as the correlation coef- 
ficient. Note that R 2 is referred to unity, thus when R 2 is 1.0, it means 
that the fitted model exactly matches the data. A low value of R 2 does 
not necessarily imply that the fitted model does not closely approximate 
the actual data, as the NIPIP technique will approximate a "best" fit of 
the data. For example, a low correlation coefficient will usually result 
when using data in which a large amount of noise is present. A low corre- 
lation coefficient should be cause for concern, and it should be 
determined whether the cause of this inability to fit the data is due to 
an inappropriate model or to problems within the input data. 


B. ADDITIONAL FEATURES OF THE NIPIP TECHNIQUE 

The previous sections described the basic function of the NIPIP esti- 
mation technique. This section presents some additional refinements or 
features of the identification procedure. 


1. Sliding Time Window 

The first of these refinements Is the use of a "sliding window." The 
sliding window Is effectively a moving average process which allows the 
ideutif ication of the time-varying behavior. Using the running solution 
for the estimation equation coefficients, the H^H and matrices will 

become "stale" and mask the time-varying behavior of the pilot when the 
curve-fitting process is to take place over a long period of time. One 
method of eliminating this "staleness" of the H^H and H^. matrices is to 
reinitialize them periodically. However, periodic reinitialization of H^H 
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and was shorn to be unacceptable In Ref. 3 because the solution for c^ 
had to reconverge after each reinitialization; thus making It difficult to 
Interpret the results. The problems of staleness and reconvergence In the 
solution can both be avoided by using the sliding window approach to up- 
dating the H*H and matrices * 

The sliding window saves N values of y and _F (the window Is thus N*T 
seconds long). After N data points have been saved and the H^H and h'jr 
computed for these N data points, the window is "slid" to include the NH 
data point and the data point at the beginning of the window Is dropped as 
follows: 


[H T H] 


n+1 


T T 

H a H - Fj F, 
n n -1—1 


+ F F 
Ai+l -n+1 


(27) 



T T 


T 

+ F \ , 
—n+1 


'n+1 


(28) 


where the vector is defined in Eq. 3. Future H T H and matrices are 
computed In a similar fashion and the process continues until the end of 
the time history. 


2. Mas Degree of Freedom 

Another feature Is the capability of adding a bias degree >.»f freedom 
(DOF) to the estimation equation. The bias DOF accounts for any uncer- 
tainties or biases In the input data, such as trimming about a non-zero 
reference operating point or unaccounted for command inputs to the over- 
all system. An example of when a bias DOF may be needed In the estimation 
equation Is given below. 

From Fig. 1, the signals which are available are the controls, _6, and 
the states, X. It is assumed 
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(29) 


6 » 


and 


& • lref-1 


(30) 


Equations 29 and 30 are used to define perturbations with respect to a 
given operating point. That is, the controls are trimmed about some oper- 
ating value, and there is some external command, X re f 9 which, in 
general, is not available for measurement. The pilot transfer function 
relates the perturbed values of the controls and states. 


6 6-6 

Y (z) » » -P -trim 

p »(*) x, X ref - X 


(31) 


cross multiplying in Eq. 31, 


^ - itrlJ D(2) ’ ( ^ref * »' 2 > 


(32) 


Moving _6_ Cr ^ in to the right-hand side and dividing both sides by D(z) yields 


6 

-P 




+ X . 
—ref 




Bias 


+ —trim 
/ 


(33) 


From Eq. 33 it is seen that the bias DOF is required when and/or 

X ref are non-zero, or their exact values are unknown. It is important to 
note that a single bias DOF can account for both of the effects of the 
unknowns X^ e f and . If J-trim an< * J^ref are known "exactly,” a bias 

DOF is not required as the differences -itrim and JCygf - JL can 

determined prior to applying the identification procedure. However, for 
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any uncertainty In the knowledge of and It is usually a good 

Idea to include the bias DOF in the estimation equation. 


3. Hiltiple Control Structures 

The AS software also has the ability to perform the least-squares 
estimation technique for several estimation equations simultaneously. The 
feature relies upon the user defined subroutine MDDE which contains the 
NAS variable specifications for each estimation equation. This subroutine 
is located within a loop in the main program. The NAS variables for each 
estimation equation are indexed within the subroutine according to the 
loop variable IM. That is, when IM * 1 the control loop specification for 
the first estimation equation is used, when IM * 2 the control loop speci- 
fication for the second estimation equation is used, and so on. An 
example of using this feature is presented in Section IV. 


4. Preprocessing of Input Data 

This feature relies upon the user-defined subroutine PREPR and is 
useful in manipulating the raw input data into the form required by the 
estimation equation. When multiple control loop structures are to be used 
and the raw input data is to be processed differently for each structure, 
the processing for each structure is again Indexed using the loop var- 
iable, IM. Thur- the same raw input data can be used for multiple control 
loop structures, where each estimation equation Is made up of different 
combinations of the same raw input data. An example of using this feature 
is presented In Section IV. 


5. Frequency Response 

When using NIPIP to identify a continuous system (l.e., a system de- 
scribed by differential equations) by means of a discrete system (l.e., a 
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z-transform of the continuous system), the identified coefficients, c, 
taken by themselves offer only limited insight into the continous sys- 
tem. Two reasons for this limited insight are that the coefficients 
change both with the data sample time and the number of degrees of freedom 
allowed for in the discrete representation of the continuous system. Thus 
it is not always possible to manipulate the identified discrete coef- 
ficients, c, back into their respective continuous system counterparts. 
In fact it is hardly ever possible to perform this reverse manipulation, 
and it is likely to be an immeuse and time-consuming effort with very 
limited results. 

This problem of Interpreting the identified coefficients can be cir- 
cumvented by interpreting the results in the frequency domain. Effective 
values of the corresponding Laplace transform frequency domain parameters 
can, in turn, be approximated from the frequency response. Besides being 
a convenient way to interpret the results, the frequency domain provides a 
means of comparing the NIPIP results with data collected by other tech- 
niques (i.e., describing function data). It also provides a means of 
computing combined system measurements (I.e., pilot-vehicle measurements). 

Prior research using NIPIP to identify continuous systems has also 
shown that interpreting the results in the frequency domain reduces the 
sensitivity of the identified solution to changes In the sample time and 

the variations in the number of degrees of freedom allowed for in the 

discrete representation. 

The stability of the frequency response (I.e., the repeatability of 
the frequency response computed at various Intervals throughout an identi- 
fication process) has been found to be closely related to the convergence 
of the identified solution (Ref. 13). That Is, a frequency response which 
is essentially invariant in amplitude and phase throughout an identifica- 
tion process generally Indicates a convergent and possibly correct 
solution. Frequency responses with large variations in amplitude or phase 
can be attributed to one or more causes. For example, the human operator 
may be a novice to the control task and may not be using a consistent 
control strategy. On the other hand, the human operator may be using a 
control strategy that Is different from that assumed by the 
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experimenter. In any case, it cannot be overemphasized that it is up to 
the experimenter, through knowledge of the control task and a thorough 
understanding of the dynamics of the controlled element, to evaluate the 
likelihood of the correctness of the identified solution regardless of the 
variations in the computed frequency responses. 
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SECTION III 

NAS VARIABLE AND SUBROUTINE DEFINITIONS 


To operate the NAS software efficiently, it is necessary for the user 
to be familiar with the program parameters which direct the flow of its 
operation (i.e., the control flags) and those used to define the problem 
at hand (i.e*, program variables). This section defines these parameters 
and describes their use in the NAS routines. The control flags are de- 
fined first, followed by the variables specifying the time history data, 
the definition of the F-vector, the definition of the transfer function 
and frequency response, and, finally, the error constraint conditions. 
The number next to the name refers to the variable number in the namelist 
subroutine • 

The NAS program parameters can be specified in two ways. The first 
method is through direct terminal input using a "namelist" subroutine 
which asks for the variable number and then the value. This method is 
called "namelist input," and the variable numbers required for this method 
are listed below. The second method is that of specifying the variables 
and their values in the user-defined subroutine, MODE. Both means of 
specifying the NAS program parameters will be demonstrated by examples in 
a later section. 

Although there are a total of fifty-seven variables listed below, most 
are defaulted *0 values that need not be changed for the majority of the 
applications. The user will, in fact, quickly learn which variables are 

relevant to a particular problem. This will usually amount to only five 
to ten variables. 

It is advisable to read and understand all of the MS parameter de- 
scriptions presented below. Whenever possible, we refer to Example A in 
Section IV when specifying the numerical value of a particular param- 
eter. After becoming familiar with the NAS parameters, the user should 
refer to the summary descriptions prese _ed in Appendix C. This will 
speed the process of setting the values of the various parameters. 
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III (Var. No. 18) 

Input control specifying the input device or logical unit 
number from which the time history data is to be read. In the 
Cyber system, this corresponds to the tape number in the pro- 
gram statement defining the filename under which the data Is 
stored. (Default: 10) 


IOOT (Var. No. 19) 

Output control specifying the output device or logical unit 
number to which the algorithm's output is directed. In the 
Cyber, this corresponds to the tape number in the program 
statement defining the filename used for output. 
(Default: 6) 


HDD (Var. No. 20) 

Made for program variable input.** 

£ 

Direct user input of NAS program variables from input 
device (e.g., terminal, etc.). 

NAS program variables are input from subroutine 
"MODE." 

(Default: 0) 


WPS (Var. No. 21) 

Specifies the number of control loop structures to be run 

simultaneously. (Default: 1, Maximum number ■ 10 on Cyber, 2 

on Apple) 

**N0TF: In order to run more than one control loop structure simul- 

taneously (i.e., NPS > 1), the program variables must be stored and 
input from subroutine "MODE;" hence MOD must be set equal to 1 for 

NPS > 1. However, a single control loop structure can be run either 
using MOD ■ 0 and NPS * 1, in which case the user inputs the program 

variables directly or using MOD • 1 and NPS * 1 and inputting the 

program variables using subroutine "MODE." 
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Direct user input is only available when running the program using 
the software under interactive processing. When using batch proces- 
sing , only subroutine MODE can be used. 


IMA (Var. No. 22) 

NIPIP analysis verification 

± 

Verification of the F-vector. This is useful to check 
that the correct F-vector has been specified. 

2 _ 

Verification of the frequency response transfer func- 
tion. This is useful to check that the proper 
transfer function has been specified. 

2 

Verifies both the F-vector and the frequency response 
transfer function. This value verifies both, whereas, 
specifying INA ■ 1 or 2 verifies only the one 

specified. 

4 


Does no verification of the F-vector or of the trans- 
fer function. This can be useful when the same 
control loop structure is being used for successive 
time histories and there is no need to verify each 
time. 

(Default: 1) 


I8A (Var. No. 23) 

System analysis (Frequency Response, Refer to Section II. B. 5, 
pp. 22-24) 

1 


Does not compute a frequency. If ISA ■ 1, the user 
can omit input of frequency response variables and set 
INA * 1 or 4 depending upon desire for verification of 
F-vector. 
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2 

Frequency response Is calculated for desired transfer 
function. Frequencies at which response Is calculated 
are output, along with amplitude ratio and phase. 

2 

Frequency response is determined for desired transfer 
function; however, a search Is made of the frequency 
response, and only those frequencies at which relative 
maxima occur In amplitude ratio are output. Subrou- 
tine "MODSH" is used for this search and its outputs 
are Mode No., frequency, amplitude ratio, phase, and 
the maximum slope in arriving at each relative maxima 
occurring within the specified frequency range. 
(Maximum of 5 relative maxima are allowed over the 
specified frequency range.) 

4 


Frequency response is determined as well as a search 
performed for relative maxima. 

(Default: 1) 


I SKIP (Var. No. 24) 

Output printing option 

2 

Full print out (useful for debugging only) 

2 

Minimum print out 
(Default: 1) 


IMFW (Var. No. 25) 

Time window options (Refer to Section II.B.l, pp. 19-22) 

2 

Sliding time window 
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ft* 

Sliding time window inoperable. 

Finite-length time window 
(Default: 1) 

**NOTE: The same result will be obtained by specifying INFW * 0 and a 

time window length greater than the length of the time history used in 
the identification as that by specifying INFW - 1. This Is useful 
when the amount of storage required for the length of the time window 
exceeds the storage allowed. 


B. TIME HISTORY VARIABLES 

These variables specify the times at which events (inputs, outputs, 
starting, stopping) or lengths of events (length of idenf iclation, length 
of sliding time window) occur within the identification process. The 
variables are usually specified in units of time, as data is usually sam- 
pled in time (e.g., samples per secor..', minute, hours, etc.); however, If 
the data is sampled in some other manner (e.g., in terms of distance, 
etc.) it is necessary only to substitute those units where the units of 
time appear in these definitions. In the following definitions it is 
assumed that the time history sample rate is in terms of samples/second. 
Thus when "seconds" appear in the following definitions we are referring 
to seconds of data on the input data file. 


T8T (Var. No. 1) 

The number of seconds of data to be ignored before starting 
the identification procedure. (Default: 0 sec) 


TEMP (Var. No. 2) 

The total length in seconds of the identification procedure, 
including the ignored time. (Default: 100 sec) 
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TOT (Var. No. 3) 

The sampling period In seconds used In the Identification 
procedure, and it may be different from the time history sam- 
ple period. However, It Is required to be In even Increments 
of the time history sample period. For example, If the tJme 
history sampling period Is x seconds and a sample period equal 
to twice this Is to be used in the Identification procedure, 
TDT should be set equal to 2x seconds . 

The sample rate used in the identification procedure should be 
based upon the expected bandwidth, u> c , of each control loop, 
with high bandwidths requiring high sampling rates ani low 
bandwldths allowing the use of low sampling rates. As a rule 
of thumb the minimum sampling rate should be: 


TDT 



[Default: 0.1 (10 samples/sec)]. 


(34) 


TOOT (Var. No. 4) 

The time in seconds between outputs of th> 1 - cification 
estimates and system analysis routines. (Dc •• . sec) 


TC (Var. No. 5) 

The time in seconds between computations of the identification 
estimates of the F-vector coefficients. This is usually equal 
to TOUT, as it is not usually necessary to compute the coef- 
ficients more (or less) often than they are output. It is 
important to note that the number of times these estimated 
coefficients are computed will not affect their values at 
future times or at the same points in time, i.e., with this 
variable equal to TOUT, it merely allows the user to monitor 
the estimation procedure as often as is desired. (Default: 5 

sec) 


TV (Var. No. 6) 

The length in seconds of the sliding or finite time window. 
This needs to be specified only if INFW is equal to zero (0) 
or two (2). The length of the time window should be based 
upon the expected bandwidth, w c , of each control loop struc- 
ture; with the higher the expected bandwidth the shorter the 
required time window and the lower the expected bandwidth the 
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larger the required time window. As a rule of thumb an in- 
itial estimate for the length of the time window Is 


TW 


10ff 

u> 

c 


(35) 


The approximate maximum length of the sliding time window is 
determined by 


^max * [(amount of storage in VAR array)/(NVAR)]*TDT 


However, if specifying a time window length greater than that 
allowed, the program will reset TW to the maximum length al- 
lowed. Hence the maximum sliding time window is easily 
obtained by specifying a time length greater than the length 
of the time history. (Default: 100 sec) (Approx, max. for 

NVAR =■ 2 and TDT « 0.1 is 12.5 sec) 


C. HIPIP VARIABLES 


This group of variables defines the estimation equation describing the 

control loop structure to be identified. These variables are probably the 
most confusing to learn for the novice user. As a teaching aide, we will 
refer to Example A presented in Section IV.A (pp. 44-61) when specifying 
values for the NIPIP variables shown below; therefore, the reader should 
review thest pages before proceeding. 


.X (Transparent to user. Refer to Section III.F, pp. 38-40) 

This is an internal array of length m, where m is the number 
of variables stored in a single record on the data input 
file. This array accomplishes the data transfer from the data 
input file to che program. 

For Example A in Section IV, A: X(l) * time, X(2) - 0, 

X(3) - > e , X(4) - w g , and X(5) - h. 


TR- 1188-1 


31 



ORIGINAL F 
OF POOR QUALi T. 

MC (Var. No. 26) 

Specifies the number of estimation coefficients on the right- 
hand side of the estimation equation (i.e., the size of the 
F-vector), that is, NC is the number of degrees of freedom. 

For Example A in Section NC * 5; aj, a 2 » b|, b 2 » 

and B. 


IBIAS (Var. No. 27) 

Specifies a bias degree of freedom (DOF) to be added to the 
F-vector. This DOF is treated as a coefficient and, as such, 
must be included when specifying NC. (Refer to 

Section II. B. 2, pp. 20-22' 

0 


Do not add bias DOF to F-vector. 


1 


Add bias DOF to F-vector. 

For Example A in Section IV. A: IBIAS = 1. 


NVAR (Var. No. 28) 

Specificies the total number of independent parameters in the 
estimation equation. An independent parameter is defined as n 
element of the _X array which is used in forming the estimation 
equation. Using an element of the X^ array more than once 
(i.e., using its current and/or past values) it is considered 
a single independent parameter. 

For Example A in Section IV .A: NVAR = 2; <S fi and 9. 


INDEX (Var. Nos. 29 to 33) 

This is an array of length NVAR which specifies the elements 
of the X^ array to be used in forming the estimation equa- 
tion. The elements of the INDEX array are specified using the 
subscript of the array element corresponding to the param- 
eters used in the estimation equation. The first element of 
this array specifies the _X array element corresponding to the 
dependent variable of the estimation equation (i.e., left-hand 
side) as well as when this array element is used as an inde- 
pendent variable in the estimation equation (i.e., when the 
F-vector contains past values of the dependent variable). The 
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succeeding elements of this * * n y specify the JC array elements 
corresponding to the independ* parameters of the estimation 
equation (i.e., right-hand side) (past values of the dependent 
or Independent variables are not considered as new variables 
in the estimation equation). The Indexes or subscripts of the 
INDEX array are given values which are used for future refer- 
ence in specifying the F-vector. The first index is given a 
value of zero (0) and each succeeding index Is given a value 1 
greater than the previous index. 

For Example A In Section IV. A: INDEX(l) = 3, X( 3) * 6 e . 
INDEX( 2) « 2, X(2) * Q 


WWtC (Var. Nos. 34 to 43) 

This Integer array of length NC* specifies how the independent 
parameters, as well as their past values, are entered into the 
F-vector. The arrangement of the parameters in the F-vector 
is determined by the order they are specified In the NVEC 
array. The elements of the NVEC array require an input of the 
form KI, where I points to the parameters in the INDEX array 
and X specifies the current or past values of the parameter. 
The codes for using I and K are shown below: 


K 


Past Value 


0 Use current value 

1 Use first past value 

2 Use second past value 


K Use Kth past value 


I (Reference value) Points to 


0 Dependent parameter, i.e., INDEX(l) 

1 First Independent parameter, i.e., INDEX(2) 


I Ith independent parameter, i.e., INDEX(I+1) 


*NC-1 if a bias DOF is used. 
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For example, NVEC(5) ■ 32 would be used to specify that the 
fifth element of the F-vector is the third (K**3) past value of 
the third parameter in the INDEX array (i.e., the second ( 1*2 ) 
Independent parameter) • 

For Example A in Section IV. A: NVEC(l) = 10, S e ; 

NVEC(2) - 20, 6 e ; NVEC(3) - 11, 9 n _ t ; NVEC(4) * 21, e n _ 2 n .~l 
n— 2 

■? (Var. Nos. 13 to 17) 

This array of length NVAR specifies the reference value as- 
sociated with each parameter defined in the INDEX array. The 
first element of BV specifies the reference value for the 
first element defined in the INDEX array and so forth. 

When using the user-defined subroutine MODE, BV is a two- 
dimensional array where the number of rows equals NVAR and the 
number of columns equals NPS . For example, BV(3,2) would be 
used to specify the reference value for the third element of 
the INDEX array defined in MODE for control loop structure 
number two. 

For Example A in Section IV. A: BV(1) * -.10681, Trim <5 e 

(rad); BV(2) = 0.0889, Trim 0 (rad). 


D. FREQUENCY RESPONSE VARIABLES 

This group of variables specifies the transfer function for which a 
frequency response is to be determined. The transfer function is made up 
of terms of the estimation equation and Includes the dependent variable, 
i.e., left-hand side, as well as estimation coefficients, c_. To obtain a 
frequency response the user must set the ISA flag (see Subsection A 
above ) • 

The range of the frequency response is based upon the bandwidth of the 
system to be identified. The frequency range should be specified such 
that it covers approximately a decade on a log plot of w, with u> c roughly 
centered within the frequency range. For our example, u) c = 2.5 rad/sec; 
hence a frequency range of w « 0.5 to 8 rad/sec should be sufficient. 


34 
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WOT (Var. No. 7) 

The minimum frequency in the frequency response. (Default: 
0.1 rad/sec) 

For Example A in Section IV. A: WMIN « 0.5. 


WHAI (Var. No. 8) 

The maximum frequency in the frequency response. (Default: 
10. rad/sec) 

For Example A in Section IV. A: WMAX * 8.0. 


MIC (Var. No. 9) 

The multiplying increment by which the frequency increases 
from w ffl £ n to w max . Successive frequencies are determined by 

w n - w n-l * w inc 


where w min < w n * w max < 36) 

(Default: 2.) 

For Example A in Section IV. A: WINC * 2. which gives fre- 

quencies of 0.5, 1, 2, 4, and 8 rad/sec. 

NOTE: The frequency response is calculated for a maximum of 25 

frequencies between w m ^ n and w max ; hence if w^ nc is too small and 
it requires more than 25 frequencies as defined by Eq. 36 to cover 
the frequency range only the first 25 frequencies will be used in 
calculating the response. 


B PC (Var. No. 10) 

A factor in the transfer function numerator which multiplies 
the high frequency gain by the value specified. Usually this 
factor is 1 or -1 , which changes the phase such that it lies 
between 0 and -180 deg. The frequency response phase angle as 
output by the NAS program will lie between 90 and -270 deg. 
Thus, in those cases where the phase line crosses these boun- 
daries, it may be difficult to interpret the phase with large 
discontinuities (e.g., going from +89 deg to -269 deg), and 
changing HFG from 1 to -1 or vice versa will change the phase 
by +180 deg, and it will then lie within the phase boundaries. 
(Default: 1.0) 
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For Example A In Section IV. A 2 riFG - 1. 
m (Var. No. 44) 

Specifies the order of the numerator of the transfer function, 
Its value is equal to the number of estimation coefficients of 
the c vector in the numerator. 

For Example A in Section IV. A: NB = 2; bj and b 7 . 


MA (Var. No. 45) 

Specifies the order of the denominator of the transfer func- 
tion; its value is equal to the number of elements of the 
estimation equation in the denominator. 

For Example A in Section IV. A: NA * 3; 1, a 2 . 


PAP (Var. No. 46) 

Specifies a pure time delay in the transfer function. Its 
value is equal to the number of samples of delay. KTAU must 
be greater than or equal to zero, as there is no prediction of 
information prior to the current value. (Default: 0) 

For Example A in Section IV. A: KTAU * 0. 


MBI (Var. Nos. 47 to 51) 

This array of length NB specifies which coefficients of the 
F-vector make up the numerator of the transfer function. The 
coefficients are usually specified in increasing powers of z 
(i.e., the lowest power of z In the numerator is the first 
element of NBI) . 

For Example A in Section IV. A: NBI(l) * 3, b,; 

NBI( 2) * 4, b 2 . 

MAI (Var. Nos. 52 to 56) 

This array of length NA specifies which elements of the esti- 
mation equation make up the denominator of the transfer 
function. When specifying a transfer function denominator 
which contains the dependent variable of the estimation equa- 
tion, i.e., the left-hand side, the first element of WAI must 
be equal to zero. 

For Example A in Section IV. A: NAI(l) * 0, 1; NAI(2) * 1, 

aj; NAI( 3) - 2, a 2 . 
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IFOPT (Var. No. 57) 

Specifies the form of output for the amplitude ratio of the 

transfer function. 

1_ 

Gives amplitude ratio squared. 

2 

Gives the linear amplitude ratio 

2 

Gives the amplitude ratio in dB. 

(Default: 3) 

For Example A in Section IV. A: IFOPT * 3. 

E. HARMING TOLERANCES 

This group of variables specifies the warning criteria which trigger 
warning messages within the NAS software. These warning messages are 
activated when the criteria are not satisfied. The warning criteria are 
specified as the maximum or minimum allowable ratios of certain variables 
within the software. Some of the reasons why these warning messages ap- 
pear are also given in order to allow the user to determine if any 
corrective action should be taken. 


TOL (Var. No. 11) 

Specifies the maximum ratio of the largest to smallest ele- 
ments along the diagonal of the H^H matrix. The purpose of 
this warning is to alert the user to the possibility of an 
ill-conditioned H^H matrix which may affect the solution for 
the estimation coefficients. An "ill-conditioned" H~H matrix 
usually results from large differences in the numerical values 
of the variables used to form the F-vector, e.g., if one var- 
iable is in units of degrees and another is in units of 
radians. This problem can usually be avoided by scaling the 
variables used to form the F-vector or using a consistent set 
of units. (Default: 1000) 


TR- llSu-l 


37 



RLD (Var. No. 12) 


ORIGINAL PAu- 13 
OF POOR QUALITY 


T 

Specifies the minimum ratio of the determinant of the H H ma- 
trix to the smallest element along the diagonal of the H T H 
matrix. The purpose of this warning is to alert the user that 
the F-vector may contain redundant DOF, which can cause a 
fragile solution for the estimation coefficients. Redundant 
DOF result in linear dependence within the H% matrix, causing 
the determinant to approach zero. This warning is used in 
conjunction with the correlation coefficient, R , to assess 
the fragility of the estimated solution. The number of DOF 
are assumed appropriate when R^ does not increase with ad- 
ditional DOF, and/or this warning does not appear. 
(Default: 0.0001) 


F. USER DEFINED SUBROUTINES 


This section defines the two user-defined subroutines which may be 
tailored to suit any particular estimation scheme. 


MDDE (involved by control flag, MOD; Refer to Section III. A, 
pp. 26-29) 

Specifies the NAS variables when several control loop struc- 
tures are being implemented simultaneously. It is necessary 
to specify only those NAS variables which are not defaulted or 
change from one control loop structure to another. To access 
this subroutine the control flag, MOD (variable No. 20), must 
be greater than zero. The default listing of MODE is given in 
Table 2. For additional examples of using subroutine MDDE see 
Section IV. 


PREPR 


Performs any processing of input data prior to use in the 
NIPIP algorithm. This subroutine transfers the data from the 
data input array, DI, to the internal array, X, which is used 
when specifying the NIPIP variables. (Refer to Section III.C, 
pp. 31-34.) The default and most common use of subroutine 
PREPR Is to simply transfer the data input file into the X 
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TABLE 2 

DEFAULT LISTING OF THE MODE SUBROUTINE 


5.5 SUBROUTINE NODE 

THIS IS A USER DEFINED SUBROUTINE FOR SPECIFYING THE 

NAS VARIABLES HHEN SEVERAL CONTROL LOOP STRUCTURES ARE TO BE 

IMPLEMENTED SIMULTANEOUSLY 

TC ACCESS THIS SUBROUTINE CONTROL FLA6 NOD NUST BE GREATER 
THAN ZERG 


THIS IS THE DEFAULT LISTING 

SYSTEMS TECHNOLOGY. INC. NAY 1982 GDH 

SUBROUTINE NODE (IN) 


IN... SPECIFIES WHICH CONTROL LOOP STRUCTURE IS BEING USED 
NAS VARIABLES ARE SPECIFIED USING THE NAMES IN THE CONNON BLOCKS 


DIMENSION NVECIIO) ,1NDEI(5).NBI (5),NA! (5).BV(5.2) 


CONNON /TINE/ TST, TEND, TDT. TOUT. TC.TN. 

NNIN.NNAI.WINC.HFG.TOL.RLD.BV 

CONNON /BAS/ IN.IOUT.HOD.NPS.1NA.ISA.ISKIP.INFN. 
NC.1BIAS.NVAR. INDEX, NVEC, 
NB.NA.KTAU.NBI.NAI.IFOPT 


THE DEFAULT CASE IS FOR CONTROL FLAG N0D=0 ... 

THERFORE THIS SUBROUTINE IS NOT REQUIRED TO SPECIFY THE CONTROL 
LOOP STURCTURES 


RETURN 

END 
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array. The default listing of PREPR is given In Table 3. 
Additional uses of subroutine PREPR are for smoothing, 
reconstructing, scaling, etc., of the raw input data prior to 
use in the NIPIP algorithm. 


G. INPUT DATA PILE FORMAT 


The input data files used by the NAS software must be given in the 
format shown in the previous examples. This format is required to allow 
the program to verify the F-vcutor and transfer function. The format of 
the input data is specified below. 


INPUT DATA FILE FORMAT 


Line 1 : M; the number of input variables per line 


4 
k 

R FI 
RDC FPS2 
NU 

LANDING PROFILE Line M+2: File identifying comments 

♦ . 10001*01 * . 2SOOE *o4 * . 8372E *00 *.1495E*02y The first two entries of column one are used to 

compute the sample rate of the data stored on 
the input file 


Lines 2 through (M+l): The variable names and units; the first 

four characters are used to verify the F-vector and transfer 
function 


*.2000t*01 *.274QE*04 *.853 jE* 00 *.1479E*0i 
♦ . 3000€ *01 * . 2720E *04 + . Bfe94t +'--‘ *.!4o3t 1 
*.4ul'0E*(*l *.269C'€+04 * . 1 191E*0i ‘.»3l2t*Gl 
*.5000E*vl * . 2650E * O* *. 12CGE *01 *.13Gvt'G. 
+.MO0E+O1 *.2600E*('4 ♦. 11?1E*01 *.1297E*02 
♦.7000E+01 *.2 j 70E*G4 *.9016E+00 ♦. 14lBE*0: 
+.8000E+01 *.2530E*04 ♦.8694E+00 ♦. 1428E*02 

♦ .9000E*01 * . 2480E *04 *. 120BE*O1 *. 1244E+02 
* . 1000E+02 * . 2450E *04 *.9014E*G0 ♦. I395E+02 
*.1100E*02 +.2410E*04 *.99B2E*00 *.I342E*02 

♦ . 1200E*02 *.2360E*04 *.1254E*(U *.1234t*02 
♦. 1300E+02 *.23206*04 ♦. I127E+01 *.1272E*02 
♦. I400E*02 *.2290E*04 *.1352E*01 *.1192E*02 


Formatted input data: The data is 

read using a list-directed read 
statement; thus the data file cannot 
be stored as a binary data. The format 
specification need only be non-integer 
as the data is read into a real 
variable array. 


-1 .0 


0.0 


0.0 


0.0 


End-of-file marker 
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As can be seen, the last line contains a negative number in the first 
column. This is a signal to the program that it has reached the end of 
the input data file and will keep it from attempting to read past the end 
of the file. Thus the first column of the input data file cannot contain 
any negative numbers prior to the last line, or the program will terminate 
reading data at the point it encounters the negative number in the first 
column . 

Also it should be noted that the first two entires of the first column 
determine the sample rate of the data stored on the input file, and this 
sample rate is utilized when computing the number of samples to be used or 
ignored when specifying the time history variables. For example, the 
sampling rate in this case is 

2-1 * 1 sample/unit 

where a "unit" is generally one second. Since specifying TST, TEND, TDT, 
TOUT, TC, and TW in the number of units (e.g., seconds) determine the 
number of units to be used or ignored. If an incremental counter or clock 
signal is available, placing it in the first column is recommended, since 
the first column cannot contain negative numbers prior to the last line. 
However, it is not necessary that this first column be the incremental 
counter or clock, only that the user recognize that the first two entires 
of column one determine the sample rate, as used by the program, that the 
time history variables should be specified based upon this sample rate, 
and that this column cannot contain negative numbers prior to the last 
line • 
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C 5.6 SUBROUTINE PREPR 

C 

C THIS IS A USER DEFINED SUBROUTINE FOR PREPROCESSING THE INPUT 
C DATA TO NIPIP 

C 

C THIS is THE DEFALT LISTING FOR TRANSFER I NO DATA FROH THE Di 
C ARRAY TQ THE I ARRAY 

C 

C SYSTEMS TECHNOLOGY. INC. HAY 1932 GDH 

C 

SUBROUTINE PREPR(DI.X.IH.H) 

C 

C 

C DI... INPUT DATA ARRAY FROH DATA FILE. TAFE, OR FLIGHT TEST 
C I ...OUTPUT DATA ARRAY AFTER PROCESSING 
C IH... SPECIFIES NHICH CONTROL LOOP IS BEING PROCESSED 

C H . . . NUHBER OF VARIABLE S PER RECORD ON THE INPUT DATA FILE 

C 

c 

c 

DIMENSION Dhl).X(l) 

C 

c 

C TRANSFER DATA FROH DI TO l ARRAY 
C 

DO 1 1=1, H 
i«n*Dim 
1 CONTINUE 
C 
C 

RETURN 

END 

A- 
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SECTION IV 

SOW ILLUSTRATIVE EEAWLES 


The four examples presented below are designed to illustrate how to 
use the NAS software package as well as to demonstrate Its capabilities. 
Each example is presented as follows: 


1. Statement of the problem and objectives 

2. Consideration of NIPIP parameters 

a. F-vector 

b. Time window and frame time 

c. Frequency response 

3. Definition of the NAS variables 

4. Description of the NAS outputs 

5. Discussion of the results 


The first example is a single-input, single-output compensatory track- 
ing task. The second example is an extension of the first, an outer loop 
has been added to the control task, thereby providing a more complex ex- 
ample. The third example identifies the pilot's control strategy during a 
flare to landing. The fourth example Identifies a helicopter pilot's 
control strategy while performing a deceleration to hover. This last 
example demonstrates how NIPIP can be used for solving nonlinear control 
problems. 



A. A SINGLE-LOOP EXAIPLK 
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1. Statement of the Problem and Objectives 

The functional block diagram shown in Fig. 4 summarizes a single- 
input, single-output, compensatory tracking control task. The "pilot," 
depicted as Y*, is regulating pitch attitude, 6, using the elevator, 6 e . 
The aircraft is flying through atmospheric turbulence, Wg, and the compen- 
satory control task is to maintain the commanded value of pitch attitude, 

v 

The pilot, aircraft, and vertical gusts were all simulated on a digi- 
tal computer. Details of the aircraft and gust simulation are presented 
tn Appendix A. The form and numerical values of the simulated pilot con- 
trol strategy, Y”, were chosen to be representative of human pilot 
dynamics measured by other techniques (Ref. 6), a..d also to give good 
closed-loop characteristics. The Laplace transform of Yp is given below, 


e 

Y (s) 

P 


6 

e 

T 

e 


- -4.87 


■j + 3.0 


rad 6 

(— 


M / A A V / > 


(37) 


a 

The Yp(s) shown in Eq. 35 yields a bandwidth of 2.5 rad/se md a phase 
margin of 35 deg. 

Figure 5 contains time histories of the closed-loop (i.e., aircraft 
and pilot) response to vertical gusts. The objective of NIPIP is to iden- 
tify the parameters of Eq. 35 using the time histories of 6 and 6 g . For 
this example 0 C is constant throughout the time history (i.e., 0 C is equal 
to the trim values of 0). 


2. Consideration of NIPIP Parameters 

a. Specifying the F-Vector . The first thing to decide is how many 
degrees of freedom are required to describe the system to be identified. 
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Figure 5* Time History Used for Example A. 
(Closed- Loop Response to Vertical Gust.) 
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Q 

For this example we know there are four parameters describing Y^: two in 
the numerator and two in the denominator, as shown in Eq. 37. Since NIPIP 
operates on discrete points in time, we express the unknown parameters of 
Yp using z-transforms (see Entry 5 in Table 1). Thus, 


8 e 

Y (*) “ / 

p ®e 


. -1 . . -2 
bjZ + b 2 z 

~ -I ^2 

1 - a^z - a^z 


(38) 


where ® e “ ® c ” ®» as shown in Fig. 4. Using the notation derived in 
Section II, Eq. 38 can be used to write a difference equation of the form. 


*i®e + a 2®e + b 1®«» 

1 e n-l L e n-2 1 e n-l 


+ b 2 e e + B 


n-2 


(39) 


[ a l a 2 b i b 2 B 1 


'n-1 

S n-2 

S n-1 

*n-2 


5 In c 


The "B" term in Eq. 39 represents a bias DOF and is included for the rea- 
sons presented in Section II. 

In general one will not know how many degrees of freedom to include in 
the F-vector. The effects of over and under specifying the number of 
degrees of freedom are presented later in this example. This topic is 
also addressed in Ref. 3. 
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b. Frame Tine and Time Window . The frame time (sample period) and 
time window can be calculated based upon the bandwidth of the system to be 
identified, <*> c , as discussed in Section II. Using Eqs. 34 and 35, with 
u* c * 2*4. 


TDT * -rs — * 0.13 sec (use 0.10 sec) 

10<a> 

c 

TW * * 12.6 sec (use 11.6 sec*) 

0 ) 

c 

In general one will not know the exact bandwidth of the system to be 
identified. However estimates of <a> c can be computed based upon a closed- 
loop analysis of the task. 

c. Frequency Response . The range of the frequency response is also 
based upon the bandwidth of the system to be identified. The frequency 
range should be specified such that it covers approximately a decade on a 
log plot of a), with c*> c roughly centered within the frequency range. For 
u> c * 2.5 rad /sec, hence a frequency range of o> = 0.5 to 8.0 rad /sec should 
be sufficient. 


3. Definition of the HAS Variables 

We are now ready to define the NAS input variables presented in 
Section III. Table 4 contains a partial list of the input file to be 
used. The time histories in this file were used to construct the plots 
shown in Fig. 4. The following subsections define only those variables 


*T„ • 11 .6 sec turned out to be a convenient number for the NAS 

software, due to the dimension of the VAR storage array. 
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TABLE 4 


f9 

QUALITY 


IMPUT FILE USED FOR EXAMPLE A 
(File name is EG1 ) 


si 

TiHE SEC t 

THE' RAD 0 

ELEV RAD 6 

85 FPS W_ 

HDOT FPS — h 

EXAMPLE I 

«.OOOO0C'*E*0O +.8890780E-01 -.1048150E+00 +.0M0000E+00 +.MOOOOGE+GO 
+.1000000E+00 ♦.6899810E-01 -.1048150E+00 -.5907230E+00 -.B8G4700E-01 
♦.2000000E+00 +.8905940E-01 1047750E+Q0 ♦ - 1 955400E +00 -.9445950E-01 

+.3000000E+00 +.8931070E-01 -.1067200E+00 1904910E+01 2672200E+00 

+.4000000E+00 ♦.B95E3JOE-01 1065800E+00 ♦. 1047440E+00 3803560E+00 

♦ . 5000000E+00 ♦.8942340E-01 -.1043890E+00 ♦.3889030E+00 ♦. U413&0E-01 
+.6000000E+00 +.8946960E-01 10A36&0E+00 1520230E+01 14455B0E+00 

♦.7000000E+00 «.8998520E-01 10&4000E+00 -.27901 ODE *01 -.4754900E+00 

♦.8000000E+00 +.9044760E-01 1062470E+00 2503280E+01 -.9568670E+00 

♦.9000000E+00 ♦.910S760E-01 105S820E +00 - . 1503B20E+01 - - B341 710E+00 

♦ . IOOOOOOE+Ol +.9080740E-01 1055110E+00 +.6070840E+00 3153000E+00 

*. UOOOOOE+Ol + . 9013000E-01 10 j4o30E*00 + . 1 13178GE+01 +.2374910E+00 

+ .1200000E+01 +.8933970E-01 1059240E+00 +.5167260E+OO ♦. 4M5?OOE+O0 

♦.I300000E+OI +.8875460E-01 -.10471ME+00 +.I644180E+00 +.3705290E+00 
♦. 14000y0E*0I ♦.8840140E-01 -.1075080E+00 +.1443370E+00 +.2430000E+00 

♦. 1500000E+01 +.8828550E-01 1080340E+00 ♦ . 1 737390E-01 +.8440780E-01 

♦. 1600000E+01 ♦.B845190E-01 1081560E+00 19&2770E+01 3492200E+00 

♦. 1700000E+01 + .8948780E-01 1077570E+00 2470100E+01 -.9154460E+00 

♦. 1800000E+01 ♦. 9003320E-01 -.1068590E+00 2185740E+00 -.7933200E+00 

♦.1900000E+01 ♦. 9022900E-0I 1059060E+00 1230040E+01 -.4952280E+00 

+.2000000E+01 ♦. 90342&0E-01 1052800E+00 9719870E+00 -.50B9950En<0 

+.2100000E+01 ♦. 9050010E-01 1050940E+00 -.2495250E+01 -.5965460E+00 

♦.2200000E«01 +.9065720E-01 1052240E+00 -.1474700E+01 -.70079BOE+00 

+.2300000E+0! ♦. 9080150E-01 1055210E+00 2005610E+01 7104950E+00 

+.2400000E+01 ♦ . 9087300E-01 1058350E+00 1057210E+01 -.7187B4UE+00 

+.250C000E+01 +.9079450E-01 1060850E+00 1484380E+01 -.5534400E+00 

+.2600000E+01 +.9055360E-01 1062740E+00 -.7274750E-01 4I38530E+00 

+.27O00C0E+01 +.8995140E-01 - 
A- 

t 0 & v h 

g g 

(sec) (rad) (rad) (fps) (fps) 
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which are not defaulted. The "direct input inode" was used to define these 
variables in the NAS program. 

a. NIPIP Variables (Refer to Section III.C, pp. 31-34). From Table 4 
we can infer the X-vector, 


X - (t, 0, « e , w g , h) 


From Eq. 39 and X, (Refer to Section III.C., pp. 31-34) 


NC 

* 

5 

Number of degrees of freedom 

I BIAS 


1 

Bias degrees of freedom 

NVAR 

m 

2 

Number of independent parameters in the 
estimate equation; $ e and 6 

INDEX 

* 

3,2 

Location of independent parameters in JC 
array; X(3) - 6 e and X(2) « 8 

NVEC 

m 

10, 20, 11, 21 

How independent parameter is entered in 
the F-vector 

BV 

as 

-.10681, 0.08890 

Reference (trim) value of the independent 
variable 


b. Time History Variables (Refer to Section III.B, pp. 29-31) 


TDT 

=* 0.10 

Sampling period 



TW 

- 11.6 

Seconds of sliding or finite 

time 

window 

TEND 

* 30.0 

Total time 



TOUT 

= 5.0 

Seconds between outputs of 
estimate 

identification 

TC 

* 5.0 

Seconds between computations 
tion estimate 

of 

identif ica- 
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• Frequency Response 

Variables (Refer to Section III.D, pp. 34-37) 

WMIN 

* 

0.5 


WMAX 

s 

8.0 

♦ U) * 0.5, 1.0, 2.0, 4.0, 8.0 rad/sec 

WINC 

3 

2.0 


NB 

3 

2 

Ordgr of numerator, number of coefficients 
of c 

NA 

m 

3 

prder of denominator, number of elements of 
c 

NBI 

3 

(3, 4} 

Specifies which coefficient is in the numer- 
ator; bj and b2 

NAI 

* 

{0, 1, 2} 

Specifies which coefficient is in the denom- 
inator; 1, aj, and a£ 

1. Control 

. Flags (Refer to Section III. A, pp. 26-29) 

INA 

= 

3 

F-vector, transfer function frequency re- 
sponse option 

ISA 

3 

2 

Frequency response option 

INFW 

3 

0 

Time window option 


e. User-Defined Subroutines 


The default subroutines for MODE and PREPR will be used for this 
example (see Section III.F, pp. 32-35). The namelist option will 
be used to set all of the necessary NAS variables. 
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4. Description of the HAS Outputs 

Table 5 contains an annotated output of the NAS program. Items typed 
by the user are underlined. The outputs for both 5 and 7 DOF are Included 
In the table. 

A 

Note that the solution for _c_ Is convergent after only 50 data points 
(see output at t * 5.4 sec), and Is Invariant throughout the run. For the 
7 DOF case, a "warning" message Is output due to an "overspecified 
F-vector." This message resulted from the determinant of the H T H matrix 
being very small (i.e., 10“^), which indicates that there are redundant 
degrees of freedom in the F-vector. 


5. Discussion of the Results 

Figure 6 contains frequency responses of the exact and estimated pilot 

0 6 
describing function, Yp(jw). Note that the estimated and exact Yp(j<*>) are 

In perfect agreement for the 5 DOF case. 

A A 

In addition to the 5 DOF case, Fig* 6 also presents the Yp(jo>) using 
3, 4, and 7 DOF. The correlation coefficients, R^, and determinants of 
the H*H matrix, A, are also shown in the figure. Note that the 3 DOF case 
did not fit the true frequency response and that is 0.88. The 4 DOF 
case caua closer to fitting the true frequency response, and increased 

to 0*98. The 7 DOF case produced the same result as the 5 DOF case, with 
* 1.00, but note that DHTH is 0.1 * 10~^! This low a value of A indi- 
cates a linearly dependent H A H matrix which makes the solution for Yp very 
’'fragile." That is, any noise In the measurements will cause variations 

f 5*. 

A Q 0 

The logical choice for Y” (even If we did not know the true Y”) would 
be the 5 DOF case. This choice Is based both on the fact that peaks at 

5 DOF and the dramatic decrease In A In going from 5 to 7 DOF. 
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NAS OUTPUT FOR EXAMPLE A 


NAS 


Execute the NAS program 


INPUT FILE 1 EG1 Filename 


1 

TIME 

SEC 

4 

THEI 

RAD 

3 

ELEV 

RAD 

4 

NS 

Fr'S 

5 

HDOT 

FFS 


FILE ID : EXAMPLE 1 NO. VAR. 5 

DIRECT INPUT USINS NAMELIST 


VAR. HO.? 26 


NC ~ 

0 CHANGE TO 5 

NC 

5 

VAR. MO. 7 n 


I3IAS = 

0 CHANGE UM 

I BIAS 

1 

VAR. NO. 1 28 


NVAR 

0 CHANGE T0_2 

NVAR 

2 

VAR. NO.’’ 29 


INDEX 1 = 

0 CHANGE TO 2 

INDEX 1 = 

J 

VAR. NO. 1 ;o 


INDEX2 = 

0 CHANGE TO 2. 

INDEX2 = 

2 

VAR. NQ.’M 


NVEC1 

0 CHANGE TO JO^ 

NVEC1 = 

10 

VAR. NO.’ 05 


NVECI = 

0 CHANGE TO 20 

NVEC2 

20 

VAF. NO. - ’ 3« 


N*EC3 = 

0 CHANGE T0_U 

N E r 3 = 

i t 

. <r. mo.' 


NVECI 

(■ CHANGE TO 21 

NVEC4 - 

21 

VAR. NO. 1 J2 


BVI 

V.voO CHANGE TO -.10681 

BV1 

-. 1063 

VAR. MO . 1 J_4 


BV2 

0.000 CHANGE TO .0089 

BV2 

.8890E-01 

VAR. NO. 1 3. 


TDT “ 

. 1000E+00 CHANGE TO .1 

TDT 

.1000 


Note: User inputs 

are underlined 


t 

Use namelist to 
set NAS variables 
for 5 DOF case 
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VAR. HO.? 7 


WHIN 

, 1000E+00 CHANGE T0J5, 

WHIN 

.5000 

VAR. N0.’.w 


WHAX 

10.000 CHANGE TOj). 

MAI 

8.000 

VAR. NO.M 


WINC 

2.000 CHANGE TOJ^ 

NINC - 

2.000 

VAR. NO.’ 44 


NB 

0 CHANGE TO 2 

N8 

2 

VAR. NO.’ 45 


NA s 

0 CHANGE T0_3_ 

NA 

*> 

J 

VAR. NO.’ 47 


NBU 

0 CHANGE TO 2 

NBI1 

3 

VAR. NO.’ 48 


NBI2 

0 CHANGE TO 4 

NBI2 

4 

VAR. NO.' 1 ^2 


NAI1 

0 CHANGE TO _0. 

NAI1 

0 

VAR. N0.’_53_ 


NAI2 

0 CHANGE TO 1 

NAI2 

1 

VAR. NO.’Ji 


NAI3 

0 CHANGE TO 2 

NAI3 

2 

VAR. NO. ’21 


INA 

1 CHANGE TO 2 

INA = 

3 

VAR. NO. ’21 


ISA 

1 CHAN6E TO 2 

ISA 

2 

VAR. N0.”_25 


INFW 

1 CHANGE TO £ 

INFH 

0 

•.’iB Mn. 1 -i 



I (DEPENDENT VARIABLE* , FOR CONTROL LOOP 1 

0 ELEVNN-O) 

THE F-VECTOR (INDEPENDENT VARIABLE): 

1 ELEV(N-l) 

2 ELEViN-2) 

3 THET(N-l) 

4 THET (N*2) 

5 BIAS DOF 


ORIGINAL PAGE IS 
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! A negative subscript 
signals the end of 
the namelist input* 



l - [». 


n-1 n-2 


n-1 


W' B ) 


* A zero subscript causes all of the namelist variables tc be output. 
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TABLE 5 (Continued) 


original 

0F POOR 


PAGE 13 

Quality 


TRANSFER FUNCTION FOR CONTROL LOOP 1 ) 5 

NUN: THET(N-l) THETIN-2) 211- 0 >— (z) 

DEN : 1.0 ELEV(N-l) ELEVIN-2) ) e 


-1 * -2 

c,z + c. z 

3 4 

~ -1 ^ -2 

T - c.jZ - c 2 z 



VERIFY VARIABLES TO BE SAVED 
ELEV WTH TRIS VALUE OF -.10681 

THET NITH TRIM VALUE QF .06890 


Note that NAS adjusted T w 
of 11.6 sec.* 


TN, TOT. TC. TOUT = .116E+02 .100E+00 .500E<01 .500E+01 


to a maximum 


END VERIFICATION... If F-vector and transfer functions 

TYPE 'O' TO CONTINUE; M* TO RETURN TO NAMELIST? ± are okay, type 0, otherwise type 

-1 and correct problem 

T 1(1) Y2B DHTH R2 N CONTROL LOOP 
5.4 5.3 . 143E-03 .9B2E-1B .100E+01 50 1 


CHAT- . 14042E+01 -.697S3E*00 .44813E+00 -.32906E+00 -.25216E-05 
FREQ / AMP / PHASE 



T 

10.4 


.50 

1.00 

2.00 

4.00 

8.00 

•7.67 

-7.20 

-5.49 

-.20 

.33 

4.87 

9.18 

14.46 

5.00 ■ 

*88.44 

X(l) 

m 

DHTH 

R2 

N 

10.3 

.227E-05 

.913E-15 

.1006 '91 

100 


oh (rad/sec) 

4 Y p (jco) 

CONTROL LOOP 
1 


IXp (* o )| 


CHAT- .14043E+01 -.69764E+00 .44812E+00 -.32905E+00 -.249696-05 
FREQ / ANP / PHASE 


.50 

1.00 

2.00 

4.00 

8.00 

-7.67 

-7.20 

-5.49 

-.20 

.33 

4.37 

9.18 

14.46 

5.00 

-38.44 


*The maximum allowable is a function of the memory available on the 
computer used by NAS. 
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TABLE 5 (Continued) 


ORIGINAL P*G£ U 
OF POOR QUALITY 


T lit) Y28 DHTH R2 N CONTROL LOOP 
15.4 15.3 . 356E-05 .797E-14 ,I00E*01 116 1 


CHAT: ,14G43E*01 -.6?768E*00 . 4480SE«-00 -.32898E*00 -.25137E-05 


FREQ / ANP * PHASE 



.50 1.00 

2.00 

4.00 

8.00 



-7.67 -7.20 

-5.49 

-.20 

.33 



4.37 9.17 

14.46 

5.00 

-88.45 


I 

111) Y2B 

DHTH 

R2 

N 

CONTROL LOOP 

20.4 

20.3 . 397E-05 

1 

U1 

oo 

.100E+0I 

116 

1 

CHAT* 

.14043£K>1 

- . 69773E+00 

. 44801E*00 -.32893E+00 -.24797E-05 

FREQ / 

ANP / PHASE 






.50 1.00 

2. JO 

4.00 

8.00 



-7.67 -7.20 

-5.49 

-.20 

.33 



4.86 9.17 

14.45 

5.00 

-88.46 


T 

1(1) 72B 

DHTH 

R2 

N 

CONTROL LOOP 

25.4 

25.3 .544E-05 

.280E-13 

, 100E+01 

116 

1 


CHAT= . I4044E»01 -.69775E+00 .44800E+00 -.32893E+00 -.24976E-05 
FFES < AflP / PHASE 


.50 

I. JO 

2. 'JO 

4.00 

8.00 

-7.6' 

-\20 

-5.49 

-.20 


4.36 


14.45 

5.00 

-33.46 


I Last print out 
for 5 DOF 
case . Note that 
solution was invariant 
throughout run 
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TABLE 5 (Continued) 


ORIGINAL P,%2- fj 
OF POOR QUALITY 


INPUT FILE ” JG1. 

1 TINE SEC 

2 THE! RAD 

3 ELEV RAD 

4 Ho FPS 

5 HDOT FPS 


Program ].oops bacic 
to this po’ . « All 
NAS variab.* axe 
unchanged from above 


FILE ID : EXAMPLE I NO. VAR. 5 


DIRECT INPUT USING NAMELIST 


VAR. 

NO. 

‘I 




TEND 


Z 

6.1 

300 

CHANGE 

TEND 


z 

12 

.00 


VAR. 

NC. 

"ii 




NC 


= 

5 

CHANGE 

T0_7_ 

NC 


= 

7 



VAR. 

NO. 





NVECI 

/ 

= 

0 

CHANGE 

T0_30_ 

NVEC5 

= 

30 



VAR. 

NO. 

’ j? 




NVEC& 

= 

0 

CHANGE 

TO _31_ 

NVEC6 

- 

31 




NO. 





NB 


- 

T 

L 

CHANGE 

702 

NB 


= 

7 

J 



VAP. 

NO. 

1 45 




NA 


= 

T 

J> 

CHANGE 

T0_4_ 

NA 


= 

4 



VAR. 

NO. 





N8I3 


= 

0 

CHANGE 

TO 2 

NB 1 3 



5 



VAR. 

NO. 

’JB 




NAM 


s 

0 

CHAN6E 

T02 

NAM 


= 

6 



VAR. 

NO. 

’ -i 





r 

Use namelist to 
set NAS variables 
for 7 DOF case 


V (DEPENDENT VARIABLE). FOR CONTROL LOOF 1 

0 ELEV IN-0 1 

THE F-VECTOR (INDEPENDENT VARIABLE): 

1 ELEVIN-li ^ 

2 ELEV(N-2) 


3 THETiN-l) I 

4 THET ( 54-21 > 

5 ELEViN-3) 

6 THETIN-3) 

7 BIAS DOF 


F = 



y = 


6 


n 


'n-2 


,6 n-r G n-2' & <= 


,0 


n- 


n-;> 
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ORIGINAL PAGE IS 

TABLE 5 (Continued) OF POOR QUALITY 


TRANSFER FUNCTION FOR CONTROL LDOF 1 

NUN : MTIN-1) THETiN-2) ELEVIN-3) Ztl- 0 
DEN : 1.0 ELEV(N-l) ELEV(N-2) THET1N-3) 


Note incorrect 
specification of 
transfer function 


VERIFY VARIABLES TO BE SAVED 
ELEV MITH TRIH'JALUE OF -.10631 

THFT NITH TRIH VALUE OF .08890 

TM.TDT.TC, TOUT * .116E+02 .i00EM)0 .500E*01 .500E+01 

END VERIFICATION... Return to r jnelist 

TYPE ’0* TO CONTINUE: M* TO RETURN TO NATEl’ET? to correct transfer 

function specification 

DIRECT INPUT USIN 6 NAMELIST 
U AR. NO . 1 49_ 

NB13 * 5 CHAN 8 E TD_ 6 _ 

NBI3 = 6 

VAR. NQ. 1 ^ 

NAM = 6 CHANGE T0_5 

NAM * 5 

VAR. NO . 1 -1 


V (DEPENDENT VARIABLE), FOR CONTROL LOOP 
0 ELEVIN-O) 

THE F-VECT0R (INDEPENDENT VARIABLE): 


ELEV(N-l) 
EL_V (N-2) 
THET(N-l) 
THET (N-2) 
ELEVIN-3) 
THETiN-3) 
BIAS DOF 


’ = [&e n-/ Be — ,9n - ,,0 n- 2,5 e- vVy B] 


'n -2 


n-3 
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TABLE 5 (Concluded) ORIGINAL fg 

OF POOR QUALITY 


TRANSFER FUNCTION FOR CONTROL LOOP 1 ) 6 e C,s -1 +c, z" 2 +c.z" 5 

NUN : THET'N-1) THEKN-2) THETIN-3) Zll- 0 > 9 “ (z) = -4?— ;;. „ — ~s — 

DEN : ..0 ELEV(N-l) ElEVtN-2) ELEVTN-3' J e 1-^z -c z -c z' 


VERIFY VARIABLES '3 BE SAVED 
ELEV KITH TRIN VALUE OF -.10631 

THET WITH TRIN VALUE OF .08390 

TN, IDT. TC, TOUT = .116E*02 -lOOE^OO .5008*01 .500E*01 

END VERIFICATION... 

TYPE 'O' TO CONTINUE: ’-l* TO RETURN TO NRNELIST n _0_ 

T XU) Y2B DHTH R2 N CGNTROL LOOP 
5.5 5.4 . 148E-05 .351E-34 .100E*01 50 1 


CHAT- .96?30E*00 -.84101E-01 .44802E+00 -.13305E+00 - 30475E*00 
-. 14390E+00 -.357 ^0E-05 

Itt WARNING lit POSSIBLY OVERSPECIFIED F-VECTOR 
BIN. DIAG. ELEN. OF "TH= .5572E-04 DETERNINANT OF HTH= .3509E-34 

TRACE OF HTH= .752 £-22 RATIO UNIT FOR DETERNINANT = .1000E-03 


Warning message 


due to small value 
of Ih^h/ 


FREO / ANP / PHASE 


.50 

LOO 

2.00 

4.00 

8.00 

1 Same aS 

-7.67 

-7.20 

-5.49 

-.20 

. 33 

/ 5 P DOF, even with 

4.87 

9.18 

14.46 

5.00 -88.44 

/ redundant DOF 
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Further discussion of the example presented above can be found In 
Ref. 3. In addition to the effects of the number of degrees of freedom. 
Ref. 3 also examines the effects of pilot remnant (which Is, in effect, 
measurement noise), length of the sliding time window, and changes In the 
parameters of Yj>. 


B. A HJLTILOOP EXAMPLE 


1. Statement of the Problem and Objectives 


The next example demonstrates the use of NIPIP to Identify multiloop 
pilot control strategies In the presence of discrete maneuvers. The mul- 
tiloop control task Is depicted In Fig. 7. The pilot's task is to 

• . 

maintain a commanded vertical velocity, h e , by feeding back h c to the 
commanded pitch attitude, 0 e . This is the "outer-loop" shown in Fig. 7. 
The pilot is also closing an "inner-loop" by feeding back 0 c to the ele- 
vator, 6 e * The inner- and outer-loop pilot control strategies are repre- 
0 h 

sented by Yp and Yp, respectively. As in Example 1 above, the aircraft is 
flying through atmospheric turbulence, Wg. The pilot, aircraft, and tur- 
bulence were simulated on a digital computer. Details of the simulation 
are contained in Appendix A. 

* 

6 h 

The simulated pilot control strategies, Yp and Y|J, were chosen to 
yield Inner- and outer-loop bandwldths of 3.0 and 0.2 rad/sec, respec- 
tively. The Laplace transforms of y£ and Y 0 are 


y£(s) 


9 c K 


0.0005 

s 


( rad /ft) 


0 0 

Y (s) - « K, » 1.675 (rad/rad) 

p e 


(40) 


(41) 
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The closed-loop response (l.e., pilot plus aircraft) to vertical gusts and 
h c is shown In Fig. 8. Hie magnitude of Wg (not shown in the figure) Is 
essentially the same as in Example 1 (see Fig. 5). From t = 0 to 
t * 15 sec the commanded vertical velocity Is zero. At t * 15 sec, h c is 
instantaneously changed to -87.5 fps. 

The objective of NIPIP is to identify Y p and Y p using the time his- 
tories for h, 0, and of Fig. 8. This example also shows how the user- 
defined subroutines PREPR and MODE are used in the NAS program. 

2. Cbnslderation of NIPIP Parameters 

a. Specifying the F-Vector . From Fig. 7 the inner and outer pilot 
control strategies are given by 


« e • *p (»C - e ) 


(42) 


e 

c 



h) 


(43) 


0 L 

where Yp and Y” are defined by Eqs. 40 and 41 above. Since NIPIP requires 
a time-domain representation of the control strategy equation, it is nec- 
essary to express Eqs. 42 and 43 as finite difference equations. This can 

be done by using the transfer function of the pilot, Y , and a table of z- 

0 h ^ 

transforms (e.g., Table 1). For the Y p and Y p given in Eqs. 40 and 41, 


f P - K e 

continuous 


z-transform 


discrete 


K 


Y — 


z-transform 

-► 


b j T z~ l 
discrete 


( 44 ) 


( 45 ) 
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Pitch Attitude, deg 


original page (9 

OF POOR QUALITY 



Figure 8. Time Histories Used for Example 2. 
(Closed- Loop Response to la step Command of -87*5 fps) 
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Substituting Eq . 44 Into Eq. 42 and Eq. 45 into Eq. 53 the discrete 
domain representation of the control laws are given by 


6 

e 



(46) 


e 

c 





h) 


(47) 


where b Q * Kq, b| - K£, and T is the frame time. 

In the case of our simplified simulation, 6 C can be obtained; however, 
in the realm of actual flight tests or pilot-in-the-loop simulations, 0 c 
is a command within the pilot's perception and is not obtainable. For 

this reason we will assume no direct 9, information and relate 9 to 9 and 

c c 

h. 

9 C can be eliminated by combining equations 46 and 47, 


6 

e 


- v + 




(h c - h) 


(48) 


multiplying both sides of Eq. 48 by (1 - z"*) gives 


(l - z -1 ) 6 e = - K 0 e(l - Z _1 )e + V^Tz” 1 (h c - h) (49) 


Using the notation derived in Section II, Eq. 49 is rewritten: 


6-6 = - K a (e - 0 J - K ft K/T(h - h ) + B 

e e , 9 n n-1 9 n v n-1 c' 

n n-1 


( 50 ) 
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where the bias term, "B" , has been added for the reasons presented In 
Section II. 

From Eq. 50, it is possible to derive the information necessary to 
define the NIPIP Inputs to the NAS program. 

b. Reference Values . Using Eq. 50, the reference or bias values, BV 
(refer to Section III.C, pp. 31-34), can be thought of in a number of 
ways. The first way is to specify BV as the trim elevator, 6 g , and pitch 

a • 

attitude, 0, values and vertical velocity, h, equal to h c * Using this 
approach the reference values would be subtracted in the NIPIP technique 
as follows: 



6 

e trim 


- K e<( 9 n - Vi> - 


9 trJ 





) 


(51) 


Comparing Eq. 50 with Eq. 51, the bias DOF, B, would account for the trim 
values of 6 g and 9 such that the B term of Eq. 50 would be equal to 


B 


'trim 


K e 9 trim 


(52) 


A second ^ay is to note that any non-zero biases in 6 e or 0 would be 

eliminated when using the differences shown in Eq, 50; thus BV could be 

• • 

specified as zero for 6 0 and 0 but still setting BV equal to h c for h. 
This specification of BV would allow the Mas DOF to be eliminated as 
Eq. 52 would be equal to zero. 

A third way is to specify zero biases for all three variables 6 g , 9, 
and h. In this case the B term of Eq. 50 would be equal to 


B - 


¥h ft c 


( 53 ) 
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In this example, the first approach for specifying BV will be used and the 
bias DOF retained in the estimation equation. 

c. Frame Time and Time Window . The minimum frame time and time win- 
dow are calculated based upon the inner and outer bandwldths of the 
system. Using Eq. 34 and 35 


TDT 


IL 


10(3) 


0.104 sec; TDT QL - io(o.2) * 1,57 


TW 


IL 


10 * 

3 


10,5 sec; TW 


OL 


10 * 

0.2 


157 sec 


where the subscripts "IL" and "0L" refer to the inner and outer loops, 
respectively. To simultaneously identify both the inner- and outer- 
control loop elements, a frame time of 0.1 sec was used since it will 
identify both control loops. A time window equal to the length of the 
time history was chosen since it was found that a time window of 157 sec 
was not required for the identified solution to converge. 

In general the frame time and time window lengths specified by Eqs. 34 
and 35 are useful as starting points when making a first pass at iden- 
tifying the control loops, with variations about these values permitted as 
deemed necessary. 

d. Frequency Response . The frequency response of the Inner loop, Y°, 
is determined by specifying the transfer function 



For w a 3.0 rad/sec a frequency range of o) * 0.5 to 8 rad/sec is speci- 
IL u 

fled. The outer-loop transfer function, Y" , can be found by taking the 

ratio 
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( 55 ) 


For u> * 0.2 rad/sec a frequency range of w ■ 0.02 to 0.8 rad/sec is 
C 0L 

specified. 


User-Defined Subroutines. The subroutine PREPR, shown in Table 6, 


is used to compute the differences, 6 n - 9 n _j and 6 g - PREPR 

defines these differences as Elements 2 and 3, respectively, in the 
X array, h is defined as Element 5 in the X array. (For a definition of 
the X array, refer to Section III.C, pp. 31-34.) 


0 V 

The forms of the control loop elements Yp and y£ are specified using 
the subroutine MODE, shown in Table 7. The following section describes 
how subroutine MODE was defined. 


3* Definition of the NAS Variables 

In specifying the MODE subroutine, only those variables which are not 

defaulted or change from one control loop structure to another need to to 

be defined. In the following subsections, the inner loop control element, 

Yp,is referred to as "control-loop 1" and the outer loop control element, 

. 

is referred to as "control-loop 2." 

A partial list of the input file to be used by NAS is shown in 
Table 8. 

a. NIPIP Variables (Refer to Section III.C, pp. 31-34). From Table 8 
and the description of subroutine PREPR we can infer the X-vector, 


X 


Wi- V‘vr “8- 
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ORIGINAL PAGE [3 
0F P0 °R QUALITY 


C 5.7 SUBROUTINE PREPR 

C 

C THIS IS AN EXAMPLE SUBROUTINE SHOWING HOW THE INPUT DATA 
C HAY BE PREPROCESSED PRIOR TO USE IN THE NIPIP TECHNIQUE 

L 

C THIS SUBROUTINE CALCULATES THE DIFFERENCES BETWEEN THE CURRENT 
C AND FIRST PAST VALUES OF THE SECOND AND THIRD ELEMENTS OF THE 

C D1 ARRAY AND OUTPUTS THEN AS THE SECOND AND THIRD ELENENTS OF 

C THE I ARRAY, ALL OTHER ELENENTS Oi- THE DI ARRAY ARE TRANSFERED 
C TO THE I ARRAY NITH NO ADDITIONAL PROCESSING 
C 

C SYSTEHS TECHNOLOGY. INC. RAY 19B2 GDH 

C 

C 

C DI... INPUT DATA ARRAY FROM DATA FILE. TAPE. OR FLIGHT TEST 
L J... OUTPUT DATA ARRAY AFTER PROCESSING 

C IN... SPECIFIES WHICH CONTROL LOOP IS BEING PROCESSED 
C N...NUNDER OF VARIABLES PER RECORD ON THE INPUT DATA FILE 

C 

r 

C 

c 

SUBROUTINE PREPR(DM.IR.N) 

C 

D I HENS ION Dim.X(l) 

C 

C 

C 60 TO 5 ON THE FIRST PASS TO PICK UP THE PAST VALUES OF THE 

C SECOND AND THIRD ELENENTS OF THE DI ARRAY 

C 

C USE IN 1 99 TO GO TO 5 ON THE FIRST PASS 
C 99 IS USED SINCE IN WILL NOT BE EQUAL TO 99 DURING FUTURE 
C PROGRAM OPERATION. 

C FOR I M*' = 2 COHPUTE DIFFERENCES 
C 

IF i IN ,EG. 99)60 TO 5 
C 

C TRANSFER DATA FROM DI TO X ARRAY 

C 

DC 1 1=1, H 
X (I)-DI (I) 

1 CONTINUE 
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c 

c 

C COHPUTE DIFFERENCES 

C 

X (2>=X (2)-TML 
X (3)-X <3>-DEL 
C 

C GP TO RETURN IF IN<2 SINCE THE SANE DIFFERENCES ARE TO BE USED 

C IN CONTROL LOOP STRUCTURES 1 AND 2 

C PICK UP NEW PAST VALUES NNEN IN-2 

C 

IF ( IN .LT. 2) BO TO 9? 

C 

5 THL S DJ 12) 

DEL-DK3) 

C 

99 RETURN 
END 
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LISTING OF SUBROUTINE MODE 


C. 5.6 SUBROUTINE NODE 

C 

C THIS IS AN EXAMPLE OF A USER DEFINED SUBROUTINE FOR SPECIFYING 

C THE NAS VARIABLES WHEN SEVERAL CONTROL LOOP STRUCTURES ARE TO 

C BE IMPLEMENTED SIMULTANEOUSLY. 

C 

C TO ACCESS THIS SUBROUTINE THE CONTROL FLAGS HOD AND NPS MUST BE 
C SET EQUAL TO "1“ AND '2" 

C RESPECTIVELY. USING THE NAMELIST DIRECT INPUT. 

C 

C THE CONTROL LOOP STRUCTURES FOR EXAMPLE 2 ARE SPECIFIED AS 
C FOLLOWS: 

C THE INNER LOOP IS C PECIFIED WHEN IM=S 

C THF OUTER LOOP IS SPECIFIED NHEN IM-2 

C 

C ONLY THOSE VALUES NHICH ARE NOT DEFAULTED OR NHICH CHANGE FROM 

C CONTROL LOOP STRUCTURE 1 TO CONTROL LOUP STRUCTURE 2 ARE SPECIFIED. 

C 

C NAS VARIABLES ARE SPECIFIED USING THE NAMES IN THE COMMON BLOCKS 

C 

C 

C SYSTEMS TECHNOLOGY. INC. MAY 1982 GDH 

C 

c 

C IM... SPECIFIES WHICH CONTROL LOOP STRUCTURE IS BEING USED 
C 

c 

c 

SUBROUTINE HODE(IM) 


DIMENSION NVEC (10) . INDEX (5) .NBI (5).NAI (5) . BV<5. 2) 
C 

c 

COMMON /TIME/ (ST. TEND. TDT, TOUT. TC.TW. 

1 NM1N.NHAX.NINC.HFG.T0L.RLD.BV 

C 

COMMON /BAS/ 1N.I0UMDD.NPS.INA.JSA.ISKIP.INFW. 

1 NC,IEIAS,NVAR. INDEX, NVEC. 

2 NB.NA.KTAU.NB1.NAI.IFCPT 
C 

C 

GO TO M.2I.IM 
C 

C NAS VARIABLES FOR CONTROL LOOP I (INNER LOOP) 

C 
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1 CONTINUE 
NC'J 
IB1ASM 
NVAR-3 
INDEX U) =3 
INDEX (2) -2 
INDEM)-! 

NVEC ( I >=01 
NVEC (2>=12 
BV(l.l)*-. 10681 
BV (2. 1 ) = .0889 
BV(3,l)s-B7.5 

C 

TST=15. 

TDTs, 1 
TEND*30. 

T0UT*3. 

TC*3. 

C 

NHIN S .5 
(WAX *8. 

NINC--2. 

C 

NB«1 
NA=l 

NBim-i 
NAIU)-0 
C 

ROOM 
NPS*2 
INA 1 ’ 

ISA-2 

c 
c 

60 TO 99 
C 
C 

C NAS VARIABLES FOR CONTROL LOOP 2 (OUTER LOOP) 
C 

2 CONTINUE 

BV( 1.2) *-.10681 
BVI2.2K0889 
BV(3,2)=-B7.5 
WUN=.02 
HNAX-. 8 
HI NC= 1 . 5 
C 

NBI < 1 ) =2 
NA I ( 1 ) = 1 
C 

c 

99 RETURN 
C 

END 
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TABU 8 

INPUT FILE USED FOR EXAMPLE B 
(File name is STEP) 


TIME SET t 

THET RmD 6 

ELEV RAD b 

MG FES w 

HD5T FF'S - h 

DISCRETE STEF -87.5 FPS 

♦.0000000E*00 ♦. 8890780E-01 -.t06815uE+0t +.OOOOOOOE+OG *.0000000 
♦.lOOOOOOEK-O ♦.6870710E-O1 -.1068150E»00 ♦. 20068 10E -03 M945720EM* 
♦.2000000E*00 ‘.881822GE-01 -.1071350E*00 ♦.7158820E-03 *.595863QE*G(. 
♦.3000000E+00 ♦. 87554 90E 01 -.I07"*>50E+00 M313370E-02 *.846825GE*o0 
♦.4000000E+00 +.8732380E-0' . I089440E+0C M502140E-02 ♦.6301350E*00 
♦ . 500C000E *(K* +.9753650E-' '92790E+00 ♦. 125716OE-02 *.2335550E*00 

«.6000000E*00 ♦.8753170E- .viMW ♦.125108GE-02 M543760E+00 
♦. 706000 0E*00 +.8739910E-O1 1083730E+OC *.1360970E-02 *.e995320E*00 

♦.800u000€*0(» + . S742340E-01 1090360E+00 M3O1680E-02 *.602l640E*» f 

*.900000'.»E*00 *.87*35o0£-0i -. 1089450E+00 ♦..259370E-02 *.6722370E*00 
♦.lOOOOOOE+OS ♦.8763130E-01 1088690E+00 M03O010E-02 +.55383?OE*00 


t 

(sec) 


-1 


(rad) (rad) 


'fps) 


h 

( fl? 3 ) 
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NC 

* 

3 

Number of degrees of freedom 

IBIAS 

- 

1 

Bias degrees of freedom 

NVAR 

m 

3 

Number of independent 
parameters in estimation 
equation 

INDEX 

m 

3, 2, 5 

NVAR specific elements in _X 
used in estimation equation, 
A0 # A6 e , h 

NVEC 

- 

01,12 

How independent parameters are 
entered in the F-vector 

BV 

« 

-.10681, 0.0889, -87.5 

Reference (trim) value of the 
parameters 


b. Time History Variables (Refer to Section III.B, pp. 29-31). 
Fcr control loops 1 and 2, 


TST = 

15. 

TEND - 

30. 

TOUT - 

3.0 

TC 

3.0 


c. Frequency Response Variables (Refer to Section III.D, pp. 34-"‘7). 
Control-Loop 1 Control-Loop 2 


WMIN * 

0.5 

WMIN - 

0.02 

WMAX * 

8.0 

WMAX - 

0.8 

WINC - 

2.0 

WINC - 

1.5 
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Control-Loop 1 


Control-Loop 2 


NB * 1 


NB - 1 


NA - 1 

See Eqs. 48 

NA - 1 

See Eq. 48 

NBI - (1> 

and 52 

NBI - {2} 

and 53 

NAI - (0) 


NAI » { 1 } 



d. Control Flags (.tefer to Section III. A, pp. 26-29). 
For control loops 1 and 2, 

MOD = 1 i These two are input 

NPS * 2 I using direct mode 

IMA - 3 

ISA = 2 


e. Description of Subroutines PREPR and MODE . Tables 6 and 7 contain 
annotated listings of the user-defined subroutines PREPR and MODE. 


4. Description of the HAS Outputs 

Table 9 contains an annotated output of the NAS program. Items typed 
by the user are underlined. 


5. Discussion of the Results 

For this case, it was not necessary to compute the frequency response 

to obtain the gains and K£. They could have just as easily been deter- 

* 

mined directly from the c's presented in Table 9. That is, From Eq. 48 it 
is recognized that Cj and 3hould equal Kq and KyK^T, respectively; thus 
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TABLE 9 . MAS OUTPUT FOR EXAMPLE B 


INPUT file *> _&TEF_ 

1 TIKE SEC 

2 THET RAD 

3 ELEV RAD 

4 HG FPS 

5 HDOT FPS 


Filename 


Note: User inputs 

are underlined 


FILE ID : DISCRETE STEF -87.5 NO. VAR. 5 


DIRECT INPUT USINS NANELIST 
VAR. NO. 7 20 

NOD = 0 CHANCE TO J_ Use Subroutine MODE 

HOD = 1 

VAR. NO. 7 21 

NPS = 1 CHANGE T0J_ — Two control loop structures 

NPS = 2 

VAR. NO. 7 -1 


Y (DEPENDENT VARIABLE). FOR CONTROL LOOP 1 Inner loop 

0 ELEV(N-O) y = B e n 

THE F-VECTOR (INDEPENDENT VARIABLE): 

1 THET(N-O) 

2 HDOT(N-l) F = [A0 , Ah B] 

3 BIAS DOF 


TRANSFER FUNCTION FOR CONTROL LOOP 1 & 

NUN : THET(N-O) IU- 0 _ * = y e 

DEN : 1.0 B e [Z) I p 


VERIFY VARIABLES TO BE SAVED 
ELEV NITH TRIH VALUE OF -.10681 

THET NITH TRIH VALUE OF .O08 7 O 

HDOT NITH TRIH VALUE OF -87.50000 

TN. TDT.TC. TOUT = .150E+02 .100E+0C ,3C0E»01 .300E+01 

Y (DEPENDENT VARIABLE). FOR CONTROL LOOP 2 Outer loop 

0 ELEV(N-O) 

THE F-VECTOR (INDEPENDENT VARIABLE): 

1 THET IN-0) ^ T f0r , 

2 HDOT(N-l) Contro1 ^ 1 

3 BIAS DOF 


TRANSFER FUNCTION FOR CONTROL LOOP 
HUH : HDOT(N-l) Ztt- 0 
DEN • -THET(N-O) 
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VERIFY VARIABLES TO BE SAVED 
ELEV WITH TRIM VALUE OF 10681 

THET WITH TR1K VALUE OF .08890 

HOOT NITH TRIM VALUE OF -87.50000 

TU.TDT.TC.TOUT = .!50E*0? . 100E+00 .300E**'l .300E+01 

NOTE: Program will not return to Namelist when MOD > O 

T ..e I/nTH F2 N CONTROL LOG- 

Ic.t !£..> . 1 1 c E - * . 73S£ ♦C , v . I00E*OI 30 1 

CHA*- . 16’»5£*01 .8S34/E-04 . 255u2E+**0 


FRE8 / A IP / PHASE 


.50 

1.00 

2.00 

4 . 0 " 

8.00 

4.47 

4.47 

4.47 

4.47 

4.:; 

0.00 

0.00 

0.00 

0.00 

0.00 


(rad/sec) 




4Y®(*d) 


J 

A 


16 . i 18.3 . 1 16E-V-: . 733E+00 .100E+01 30 

CHAT; . Ift725c+'.-: . 33347E-04 . 25502E+0O 


f'rlS 


PHASE 


CD^ ( rad/sec )- 


.02 .03 

.04 

.07 

.10 

.15 

n - 

. 34 

c 1 
. J J 

-85.49 -85. 4S 

-85.49 

-85.49 -35.49 

-85.49 

-85.49 

-85. 4- 

- 85.49 -85 

.00 0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

. 00 . 0 . 


4 y£(Jcd) A 


I Jfib Y2B DHTH R2 N 
21.4 21.3 . 1I6E-U1 . 440E+G2 .100E+01 60 


CONTROL LOOP 
1 


Repeats above 
■TOUT intervals 


CHAT; 

• 

16757E*0l 

.851 

76E-04 

.25560E»00 

r R£e * 

AMT / 

PHASE 





.50 

1.00 

2.00 

4.00 

8.00 


4.43 

4.48 

4.40 

4.48 

4.48 


0.00 

0.UG 

0.00 

0.00 

0.00 

21.4 

21.3 

. 1I6E-0I 

.44.; 

.ft +02 . 

lOOEtOl 60 

l HAT = 

, 

16?57E*0l 

.8517 

'6E-04 , 

. 2S560E *00 


FREt / AMP / PHASE 

.02 .03 .04 .07 . 10 .15 .23 .34 .51 .V 

-85.38 -85.83 -85.88 -85.88 -85.88 -85.88 -85.B8 -85.86 -35.00 -85.90 

u . 00 0.00 -.00 0.00 0.00 0.00 .00 0.00 .00 0 .'*o 
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Control 
Loop 1 
(inner 
loop) 


Control 
Loop 2 
(outer 
loop) 


in 



TABLE 9 (Concluded) 


r 




NOTES: 


ORIGINAL PAG" f.5 
0F POOR QUALITY 


T XU) Y2B DHTH R2 N COHIROL LOOF 
24.4 24.3 . 1 14E-01 .360E+03 -100E+01 90 1 

CHAT- .14748E+0! .83948E-04 .25579E+00 


FREC i AHF / PHASE 



.50 

1.00 

2.00 

4.00 

8.00 







4.49 

4.49 

4.49 

4.49 

4.49 







0.00 

0.00 

0.00 

0.00 

0.00 






24.4 

24.3 

.114E-01 

. 340E+03 . 

I00E+01 

90 

2 




CHAT= 

•14748£*01 

.83948E-04 . 

25579E+00 






FREQ i 

t AHF* / 

PHASE 










.02 

.03 

.04 

.07 

.10 

.15 

.23 

.34 

.51 

.77 


<86.01 

-84.01 - 

86.01 

-84.01 

-84.01 

-86.01 

-84.01 - 

84.01 

-B4.01 

-84.01 


.00 

0.00 

0.00 

0.00 

-.00 

0.00 

0.00 

-.00 

0.00 

0.00 

T 

U\) 

Y26 

DHTH 

R2 

H CONTROL LOOP 





27.4 27.3 . 114E-01 .114E+04 . I00E+01 120 1 

CHAT= . 14B20EKH .83567E-04 .25429E+00 

FREQ / AHF / PHASE 


.50 

1.00 

2.00 

4.00 

8.00 

4.52 

4.52 

4.52 

4.52 

4.52 

0.00 

0.00 

0.00 

0.00 

0.00 


27.4 27.3 . 114E-01 .114E+04 . 100E«01 120 2 

CHAT- . 16820E+01 .83587E-04 .25429E+00 


FRE6 / AHF / PHASE 


.02 .03 

-84.07 -84.07 
.00 0.00 


.04 

-84.07 

0.00 


.07 

-86.07 

0.00 


.10 

-86.07 

0.00 


.15 

-86.07 

0.00 


.23 

-86.07 

.00 


.34 

-E6.07 

-.00 


.51 

-86.07 

0.00 


.77 

-86.07 

0.00 


1 ) That solution for both control loops was essenti ally 
invariant throughout run 

2) Program will loop back to the same point ns in Example A 
however, in this case no further computations were made. 
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A 

these gains could be determined by direct inspection of the c coef- 
ficients* However* usually it Is not possible to directly relate the 
coefficients, c, to the pilot control parameters* In thes cases* a fre- 
quency response from which specific piloting parameters (e*g., gain and 
phase) can be deduced will aid in recognizing the identified pilot control 
elements* 


Table 10 gives the exact and 3 DOF estimated pilot gains Kq and K£* 
These gains were obtained from the frequency response amplitude ratios and 
c coefficients of control loops 1 and 2 listed in Table 9. The frequency 
response for (control loop 2) can only Identify the gain K£T as the 
knowledge that Y^ is an integrator was used in deriving the 3 DOF estima- 
tion equation* To account for the sample frame time* T* which appears in 
Eq. 50, 20 dB was added to the amplitude ratio of control loop 2 to obtain 
K£* Note that the 3 DOF estimated and exact pilot gains are in perfect 
agreemf.it* Also presented in Table 10 is the value of the correlation 
coefficient, R , which is equal to 1* 


The use of a frequency response to deduce the pilot control parameters 

0 

is further shown using 5 DOF to determine the pilot control elements Y n 


and Y . 5 DOF were obtained by not computing t'*e differences - 6 

** n 

and 9 n - 0 n _j prior to inputing them into NIFTP, but rather and 9 n _j 


'n-1 


are added as additional DOF. The estimation equation became 


n-1 


6 


e 

n 


c, 6 + c_ 9 4 c 9 + c. h + B 

1 e , 2 n 3 n-1 A n-1 

n-1 


(56) 


where c should equal 



A 

'1 


1 

A 

'2 


K 9 

C * 

b 

m 

1 

* 

CD 


C 4 


K e K h T 


_ B _ 


_ B 


( 57 ) 
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TABLE 10 

EXACT AND ESTIMATED PILOT GAINS USING THREE DEGREES OF FREEDOM 




K e 

Exact 

(rad /rad) 
Estimated* 

K£ (rad /ft) 
Exact Estimated 

R 2 

Frequency 

Amplitude 

Response 

Ratio 

1.675 

1.675 

0.0005 

0.0005 

1 

c Coefficients 

1.675 

1.675 

0.0005 

0.0005 

1 


♦Using estimate at t * 21.4 sec in Table 9* 


TABLE 11 

EXACT AND ESTIMATED PILOT GAINS USING FIVE DEGREES OF FREEDOM 



Kg (rad /rad) 
Exact Estimated* 

(rad/ft) 

Exact Estimated* 

Frequency Response 
Amp 1 .tude Ratio 

! , 675 

f 

1.66 

0.0005 

0.00041 

A 

c Coefficients 

1.675 

1.64 

0.0005 

0.00043 


♦Using estimate at t ■ 21.4 sec, NAS output not shown. 
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Values of c listed in Eq. 57 weie determined by inspection of Eq. 50. 

0 v 

In this case the apriori knowledge of Y” and YjJ Is not used explicitly 
in deriving the estimation equation. The assumed form of Eq. 56 has 
enough DOF to identify Y® and y£ as entries 1, 2, 3, and 6 of Table 1. 

Figure £ presents the exact and 5 DOF estimated pilot describing func- 

h 0 

tions of Yp. As before the frequency response Identified Yp as a pure 

gain of magnitude essentially equal to that for # 3 DOF, The frequency 
response of yJ? on the other hand now £hows that Y^ is in fact an integra- 
tor. Although the phase angle of Yp at the lower frequencies does not 
exactly Match that of an integrator, this is due to an insufficient time 
window length. (Refer to the Section II discussion on time window length, 
pp. 19-20.) The phase angle at the higher frequencies does indeed match 

that of an integrator. Further increases in # the number of DOF do not 

0 h 

change the identified solutions for Y p and Y£, and the same effects of 

increasing the number of DOF as in Example 1 were noted. 

The coefficients, c, for 5 DOF case again can be directly related to 

« 

the pilot gains Kq and using Fq. 57. Table 11 gives the gains Kq and 
• ^ 

for the 5 DOF case as determined from the c_ vector. 

C. CLOSED-LOOP PILOT-VEHICLE EXAMPLE 

1. Statement of the Problem and Objectives 

The previous examples were concerned with identifying the pilot con- 
trol elements and ignored he characteristics of the aircraft. This 
example combines the aircraft characteristics with those of the pilot, and 
the overall pilot-aircraft characteristics are identified for the landing 
maneuver (flare) of a DC-10 aircraft. 

The pilot-aircraft characteristics for the landing flare are assumed 
to be described by a second-order model with a characteristic equation of 
the form 


h + 2£u)h + 


ui^h = 


(38) 
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where h is altitude and (*) refers to time differentiation. It is further 
assumed that this characteristic equation describing the flare corresponds 
to the response from an initial altitude, h, and sink rate, h, with re- 
spect to the terminal conditions.. Thus the flare is regarded as an 
unforced response from a set of initial conditions to a set of desired 
conditions at touchdown. 

The pilot-aircraft natura - frequency, w, damping, C, of Eq. 58 

were identified from actual flight data of a DC- 10 landing maneuver. A 
detailed discussion of the DC-10 flare and pilot-aircraft model is found 
In Ref. 14. 


2. Consideration of NIPIP Parameters 

a. Specifying the F-Vecto r. Using Eq. 58, the estimated of and C 
are determined by shifting the terms containing u) and £ to the right-hand 
side. Thus the estimation equation is given by: 


• 

-2£wh 

- w^h + 

n 

n 

A • 

A 

c, h + 

c„h + B 

1 n 

2 n 

F c 



(59) 


The ,f B" term in Eq . 59 represents a bias DOF and is included for the rea- 
sons presented in Section II. 

From Eq . 59 it is seen that the pilot-aircraft parameters cj and t can 
be determined directly from the estimated coefficients, c, that is, 


U) 



(60) 
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and 



(61) 


b. Frame Time and Time Window . The frame time and time window were 
determined by the sample period and length of the landing flare time his- 
tory giving a sample period of 0.05 eec and a time window of 6.55 sec. 
This is one case where the rules of thumb are not followed when specifying 
the frame time and time window but using what is available. 

c. Frequency Response . The pilot-"trcraf t characteristics are deter- 
mined from Eqs. 60 and 61; thus a frequency response is not required and 
is therefore omitted. 

d. User-Defined Subroutines . The default for subroutines MODE and 
PREPR are used, and the NAS variables are input using the "direct input 
mode." 


3. Definition of the NAS Variables 

Table 12 contains a partial listing of the input file to be used. The 
"direct input mode" was used to define those variables which are not de- 
faulted . 

a. NIPIP Variables . (Refer to Section III.C, pp. 31-3A.) From 
Table 12 can be inferred as 


X 


{t, h, h, h) 
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TABLE 12 

INPUT FILE USED FOR EXAMPLE C 
(File name is FLARE) 


4 

TINE SEC t 

H FT ; h 

HOOT FPS h 

HDDTFPS2 fi 

DC- 10 FLARE 

+.5200E+01 ♦. 1005E+03 -.1453E+02 -.1737E+01 
+- 5250E+01 +. 1003E+03 -.1408E+02 -.2059E+01 
♦.5300E+01 ♦.9872E+02 -.1553E*02 1673E+01 

5350E+01 +.9761E+02 -.1553E+02 1995E+01 

+ .540CE+01 + 9682E+02 i;.&6E+02 -.2156E+01 

+.5450E+01 *.9619E*02 -.1564F+02 -.2091E+01 
+ .5500E'--; *.9<76E+G2 -.le41E*u2 --2349E+01 
+.5550tnil *.94I3E+02 -.1435E*02 -.215oE*01 
+ . 5600E + i +.9346E+02 1612E+02 - 19A3E+01 

+.5450E+01 +.9287E+02 -.1426E+02 -.2670EM 
+.j70GE*01 +.9140E+02 -.1485E+02 -.2542E+01 
+ . 5750E+0 1 +.9049E+02 1724E+02 -.2059E+01 

*.5800E*01 +.8954E+02 -.1746E*02 -.2284E+01 
♦ .5850EHH ♦. 8954E+02 1&73E+02 -.2509E+01 

+.590>‘>E+01 ♦.8812E+02 -.174»t*02 --.2509E 01 
+.5950EKU +.8470E+02 18HE+ f ' ' -.2316E+01 

♦.6000E+01 *.B654E*02 17S3E+02 -.2316E*01 

+.6050E+01 ♦.8622E+02 -.1709E+02 -.1963E+01 
+.6100E+01 +.8543E+02 1715E+02 -.2091E+01 

+.4150E+01 +.8464E+02 -.1719E+02 -.1769E+01 
A> 

t h A h 

(sec) (ft) (fp-0 (fps 2 ) 
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From Eq. 59 and J? 

NC - 3 

I BIAS - 1 

NVAR - 3 

INDEX - 4, 3, 2 
NVEC » 1,2 


b. Time History Variables (Refer to Section III.B, pp. 29—31 ^ 


TDT 

- 0.05 

TEND 

- 6.55 

TOUT 

■ 1 . 

TC 

- 1. 


c. Frequency Resp o nse . This section is omitted, because it is un- 
necessary. 


d. Control Flags . Use defa.' fc settir- i. 


4. Description of the NAS Outputs 

Table 13 contains ar annotated output of the NAS program. Items typed 
by the user are underlined. 


5. Discussion of the Results 

From Table 13 it is seer, that the first t..ree estimates of c(l) are 
positive in st,n; thus computing C using Eq. 61 -■ ields a negative (un- 
stable) damping value. However, the last three estimates of c(l) are 
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TABLS 13. NAS Output for Example C 

execute the NAS program 

input FILE ” FLARE Filename 


1 TINE SEC 
I H FT 

3 HBOT FPS 

4 HOOT FPS2 

FILE ID : DC-10 FLARE 
DIRECT INPUT USING NAMELIST 


NC. VAR. 4 


VAR. N0."2_ 


TEND 

100.00 CHAN6E 

TEND 

6.550 

var. v&rz. 


TC'T 

. lOOOE+OO CHANGE 

TDT 

.5000E-01 

var. mr * 


TOUT 

5.000 CHANGE 

TOUT 

1.0000 

VAR. NO." 5 


TC 

5.000 CHANSE 

TC 

1.0000 

VAR. NO.' 1 2h 


NC 

0 CHANGE TO 3 

NC 

3 

VAR. NQ.V27_ 


1 BIAS = 

0 CHANGE TO I 

IBIAS = 

i 

VAR. NO. -1 J!8 


NVAR 

0 CKAN6E TO _3_ 

NVAR * 

J 

VAR. NO.’ 29 


INDEX1 = 

0 CHANGE TO 4 

INDEX1 = 

4 

VAR. NO.’^O 


1NDE12 = 

0 CHANGE TO ^ 

index: = 

*■; 

V^R. NO.’Jl. 


INDEU = 

0 CHAF.SE TQJL 

INDEX3 * 

2 

VAR. N0.V34 


NVECi 

0 CHANGE TO 01 

NVEC1 

1 

VAR. NO.? 35. 


N’EC2 * 

0 CHANGE TO 02 

NVECi = 

2 

VAR. NC." -i 



Use namelist to set 
NAS variables 


Use negative subscript to 
exit Namelist 
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Y (DEPENDENT VARIABLE), FOR CONTROL LOOP 1 

0 H00HN-0! y = K 

THE F-VECTOR (INDEPENDENT VARIABLE)} 

1 HDOT(N-O) 

2 H (N-0) P = 

3 BIAS DOF 



VERIFY VARIABLES TO BE SAVED 
HBDT NITH TRIM VALUE OF 0.00000 

HDOT NITH TRIM VALUE OF 0.00000 

H NITH TRIM VALUE OF 0.00000 

TN.TDT.TC.TOUT - .655E+01 .500E-01 ,100E*01 .100E+01 

END VERIFICATION... 

TYPE *0’ TO CONTINUE: M* TO RETURN TO NANELISP 

T XU) Y2B DHTH R2 N CONTROL LOOP 

1.1 6.2 .465E+CI . 249E*05 .9B6E*00 20 1 

CHAT- .4U54E+00 -.66509E-01 .10778E+G2 

T XU; >26 DHTH R2 N CONTROL LOOP 

2.1 7.2 .335E+01 .107E+07 .930E+00 40 I 

CHAT- . 3820oE*00 -.67B83E-01 .10402E+02 

T XU) 12B DHTH R2 N CONTROL LOOP 

3.1 8.2 .229E+01 .150E+QB .969c*00 60 1 

CHAT- .36911E+00 -.72305E-01 . 20575E+02 

T XU) ¥26 DHTH R2 N CONTROL LOOP 

4.1 o.2 . 175E+01 .912E+08 -897E+00 80 1 

CHAT- -.48389E-01 -.52732E ?1 .18B24E+01 

T XU) ¥26 DHTH R2 N CONTROL LOOP 

5.1 10.2 . 162E+01 .347E+09 .890E<-00 100 1 

CHAT- -.68103E-01 -.529I6E-01 .15655£*0l 

T X(l) ¥26 DHTH R2 N CONTROL LOOP 

6.1 11.2 . 144E+01 .11BEU0 ,.;1E*00 120 1 

CHAT- -.21262E+00 -.54116E-01 -.7753?E*00 Use this estimate of for computing 

o> and £ 

INFLT Fill ' Program loops back to this point. No further 

runs made. 
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negative In sign which yields a positive (stable) damping value. The 

A 

estimate of c(2) on the other hand gives a real value for the natural 
frequency for all time. 

Using c(l) and c(2) at T ■ 6.1, the estimated natural frequency and 

damping are 0.23 rad/sec and 0.46, respectively. 

O 

The correlation coefficient, R , at T * 6.1 is equal to 0.86 which 
Initially leads one to the conclusion that the estimated solution does not 

match the actual data. To determine the degree to which the estimation 

solution does or does not match the actual data, the Initial condition 
response of Eq. 58 with w and C equal to their estimated values was com- 
puted. The initial conditions from the actual time history of 

h * 100.46 ft and h « -14.53 fps were used In computing the response. 

Figure 10 presents the h to h phase plane trajectories of the Initial 
condition response of Eq. 58 with the estimated o» and C as well as the 
actual data. From Fig. 10 It appears that the estimated solution does fit 
the actual data with a reasonable degree of accuracy. Thus this low value 
of R does not appear to be due to an improper modeling approach, but 
rather due to the noise which is present In the Input data. 


0. DECELERATION TO ROVER 

1. Statement of the Problem and Objectives 

The next example Is based on the research presented in Ref. 15. The 
theoretical background and derivation of the control laws are omitted In 
the following text. The purpose here is to demonstrate how the NIPIP 
technique can be used to identify the parameters of nonlinear equations. 
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Figure 10. Phase Plane Trajectories of Landing Flare 
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Hie block diagram shown In Fig. 11 summarizes the control task. Ref- 

*• 

erence 15 derives the following relationship between deceleration, R, and 
range, R, 


R 


k 2 R 


(1 + R/A)- 


(62) 


The objective of NIPIP is to identify the parameters k and A based on a 
deceleration profile of R and R. 

equation 62 shows that k and A are nonlinearly related to R and R. 
Their relationship can be made linear by algebraic manipulation of Eq. 62, 


(R/R) 1 ^ 3 - c + dR 


(63) 


It is easily shown that c and d from Eq. 63 are related to k and A in 
Eq. 62 as follows: 


k - 



(64) 


A = c/d 


(65) 


The method for using NIPIP is now clear: Identify c and d in Eq. 63 by 

" 1/3 

using the pre-processor subroutine to calculate (R/R) , then calculate k 

and A using Eqs. 64 and 65. The details of how to set up the NIPIP input 
are presented below. 
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where 

k = The effective gain for the combined pilot-vehicle system 

R = True range 

Rp = Perceived range 

A = A measure of the effective size of the objec' being viewed 
R = Range rate 

= -kR/(l + R/A) 5 
R = Range deceleration 

= k 2 R/(l + R/A) 5 

Figure 1 1 . Block Diagram for Deceleration to Hover Example 
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2. Consideration of MIPIP Paraaeters 

a. Specifying the F-Vector . The F-vector is taken directly from 
Eq. 63, 


y = (R/R) 1/3 « {R,l } { ^ } 


( 66 ) 


F c 


Thus only 2 DOF are required with one being a bias DOF (i.e., the c term 
is a bias DOF). 

b. Frame Time and Time Window . For this example the data is in the 
form of ordered pairs of R and R. Time information was not recorded; thus 
the "frame time" is in terms of the incremental number shown in 
Table 14.. All data will be processed using the infinite time window 
option. 

c. Frequency Response . No frequency response calculations are re- 
quired, because k and A will be calculated from Eqs. 64 and 65. 


3. Definition of the NAS Variables 

a. NIPIP Variables (Refer to Section III.C, pp. 31-34) The X-vector 
can be inferred from the partial list of the input data file shown in 
Table 14, 


X - {n, R, R, NU} 
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TABLE llj.. PARTIAL LIST OF THE INPUT FILE 
USED FOR EXAMPLE D 


4 

N - - ■■■■-' n 

r FT R 

RDD FPS2 K 

nu ■■ not used 

LANDING PROFILE 

♦. lOOOE+Oi ♦.2800E+04 ♦.8372E+00 M495E+02 
♦ . 2000E+-01 +.2760E+04 +.8533E+00 U479E+02 
♦.30OOE+O1 +.2720E+04 +.B494E+00 ♦.J463E+02 
*.4000E*01 + . 2&90E+04 M191E+01 +.1312E+02 
*.5000E*01 +.265OE+04 M2C8E*01 M300E*02 
+ . 6000E+0 1 +.2G00E+04 ♦. I191E+01 ♦. 1297E+02 
♦.7000E+01 ♦. 2570E+04 ♦.9016E*00 U418E+02 
♦.8000E+0! ♦. 2530E+04 +.8694E+00 +.142BE+02 
♦.9000E+01 ♦. 2480E+04 ♦.1288E*01 1244E+02 

4. 1000E+02 +.2450E+04 ♦ ^Ol&E+OO M395E+02 
♦.1100E402 +.2410E+04 +.9982E*00 *. 1342E+02 
M200E+02 2360E+04 +.1256E+01 +.1234E+02 

M300E+02 ♦. 2320E+04 M127E*01 +.1272E+02 
♦. 1400E+02 4. 2290E*04 +.1352E+01 +.1192E+02 


1.0 0.0 0.0 0.0 • End- of- file marker 
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The variable "n" is the incremental number used to number the data 
points. It is not time, but note that it is a monotonically increasing 
variable, as required by the NAS program. The variable "NU" is not used 
for this example. 

From Eq. 63 and X, 


NC - 2 

IBIAS - 1 

NVAR - 2 

INDEX - 3, 2 

NVEC - 01 


b. Time History Variables . (Refer to Section III.B, pp. 29-31) 


TDT "1.0 


With the time history variables set as shown above, NAS will process every 

A 

data point f compute, and output c every fifth data point* 

c. Frequency Response Variables # (Refer to Section III.B, pp. 29- 
31) The frequency response option will not be used for this example. 

d. Control Flags . (Refer to Section III. A, pp. 26-29) Default 
values for all of the control flags will be used. 

e. User-Defined Subroutines . (Refer to Section III.F, pp. 38-40) 
The pre-processing subroutine, PREPR, will be used to calculate the 
dependent variable in Eq. 66 (i.e., the "y" term). A listing of the 
resulting PREPR subroutine is contained In Table 15. 
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LISTING OF THE PREFR SUBROUTINE TEED FOR EXAMPLE D 


C 5.6 SUBROUTINE PREPR 

C 

C THIS IS A USER DEFINED SUBROUTINE FOR PREPROCESSING THE INPUT 
C DATA TO NIPJP 
C 

C THIS IS THE DEFAULT LISTING FOR TRANSFER 1N6 DATA FROH THE D1 
C ARRAY TO THE I ARRAY 
C 

C SYSTEMS TECHNOLOGY. INC. NAY 1982 GDH 

C 

SUBROUTINE PREPRIDM.IM) 

C 

C sassscsssssssssssassassssssssssssssssssssssssssssssssssssssss: 

c 

c DI... INPUT DATA ARRAY FROM DATA FILE. TAPE. OR FLIGHT TEST 
C I ...OUTPUT DATA ARRAY AFTER PROCESSING 

C IH... SPECIFIES NHICH CONTROL LOOP IS BEING PROCESSED 

C H ...NUHBER OF VARIABLE S PER RECORD ON THE INPUT DATA FILE 

C 

C 

C 

DIMENSION Dlil).X(l) 

C 

c 

c 

C TRANSFER DATA FROH 01 TO X ARRAY 

C 

DO I 1*1.11 
nii=Diu> 

1 CONTINUE 

«=X(2) 

RDD=X<3) 

U:j*(R/RDD)tt(l./3.> 

C 

c 

c 

c 

c 

99 RETURN 

END 

C 
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4* Description of the NAS Outputs 

Table 16 contains an annotated output of the NAS program. Items typed 
by the user are underlined. The values of k and A are calculated for a 

A A A 

few values of the Recall that k and A are related to Cj = d and C 2 L 
via Eqs. 64 and 65* 


5* Discussion of the Results 

Plots of the actual data and the NIPIP curve fit are presented in 
Fig. 12. The curve was constructed by using Eq. 62 and the indicated 

values of k and A. Note that the fit is best from R * 2800 ft down to 
about R ■ 300 ft. This is because there are few data points for 

R < 300 ft. The fit in the k < 300 ft region could be improved by not 
including the data points for R > 300 ft or using the sliding time window 
approach (i.e* t use a time window that would give a AR ± 300 to 500 ft). 
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TABLE 16. NAS OUTPUT FOR EXAMPLE D 


A>BsE6p Execute the customized, version of NAS for this exanple (i.e., 

EGD has the special version of the PREPR 
INPUT FILE 01 LAND subroutine in it). 

1 N 

2 R FT 

3 ROD FPS2 

4 NU ORIGINAL PAGE IS 

FILE ID : LANDING PROFILE NO. VAR. 4 POOR QUALITY 


DIRECT INPUT USING NAHELIST 


VAR. NO. 

? 26 





NC 

s 

0 

CHANGE 

TO 

n 

NC 

= 

n 

4 




VAR. NO. 

'll 





IBIAS 

z 

0 

CHANGE 

TO 

l 

I BIAS 

- 

1 



VAR. NC. 






NVAR 

= 

c 

CHANGE 

TO 

2_ 

NVAR 

= 

n 

4 




VAR. NO. 

?_29 





INDEX 1 

s 

0 

CHAN6E 

TO 

_3_ 

INDEX 1 

- 

3 




VAR. NO. 

? 30 





INDEX2 

= 

0 

CHAN6E 

TO 

2_ 

1NDEX2 

r 

2 




VAR. NO. 

"il 





NVEC1 

= 

0 

CHAN6E 

TO 

01 

NVEC1 

= 

1 




VAR. NO. 







( s / S )]/ 3 


1 (DEPENDENT VARIABLE) . FOR CONTROL LOOP 

0 RDD l N-0) 

THE F-VECTOR (INDEPENDENT VARIABLE): 

1 R (N-0) 

2 BIAS DOF 


y = 


l - [R n ,B] 


VERIFY VARIABLES TO BE SAVED 

RDD NITH TRIM VALUE OF 0.00000 

R NITH TRIN VALUE OF 0.00000 

TN.TDT.TC.TOUT * . 100E+03 .100E+01 .500E+C1 .500E+01 

END VERIFICATION... 

TYPE *0* TO CONTINUE: ’-l' TO RETURN TO NAMELIST? 0 
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T 

XU) Y2B DHTH 

R2 

N 

CONTROL LOOP 

7.0 

7.0 . 185E+03 .7W»E+05 

•998EI00 

5 

1 

CHAT* 

.24318E-02 . 7I437E+01 




T 

XU) Y2B DHTH 

R2 

N 

CONTROL LOOP 

12.0 

12.0 . 1Q1E+03 . 132E+07 

.997E+00 

10 

1 

CHAT* 

. 24003E-02 .73195E+01 




T 

XU) Y2B DHTH 

R2 

N 

CONTROL LOOP 

17.0 

17.0 , 173E+03 .464EI07 

.998EHRI 

IS 

l 

CHAT* 

.29171E-02 .59754E+01 




T 

Mi Y2B DHTH 

R2 

N 

CONTROL LOOP 

22.0 

22.0 . liGE+03 .210E+0B 

.997E+00 

20 

1 

n 

£ 
— » 
11 

.4I334E-02 . 26491E+01 




T 

X> 1/ Y2B DHTH 

R2 

H 

CONTROL LOOP 

27.0 

27.0 . 151E+03 .501E+08 

.997E+00 

nc 

LJ 

1 

CHAT* 

. 39948E-02 .323G3E+01 




T 

XU) Y2B DHTH 

R2 

N 

CONTROL LOOP 

32.0 

32.0 . 142E+03 . 105E+09 

. 997E+00 

30 

1 

CHAT* 

. 40853E-02 . 30704E+01 




T 

XU) )2B DHTH 

R2 

N 

CONTROL LOOP 

37.0 

37.0 . 133E+G3 .195E+09 

,997E*00 

35 

1 

r, 'J- 

.40428E-C2 .31226E+01 




T 

XU) V2B DHTH 

R2 

N 

CONTROL LOOP 

42.0 

42.0 . I25E+03 . 337E+09 

.998E+00 

40 

1 

CHAT* 

.40359E-02 .31376E+01 




T 

XU) V2B DHTH 

R2 

H 

CONTROL LOOP 

47.0 

47.0 . 117E+03 .541E+09 

.998E+00 

45 

1 

CHAT* 

.40502E-02 .31073E+01 





d = 0.002432 
c = 7.144 
k = 0.0524 
A = 293*. 


d = 0.004350 
c = 3.107 
k = 0.1826 

A = 767 . 
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TABI£ 16 (Concluded) 


T 

t<n 

<2B 

DHTH 

R2 

N 

CONTROL LOOP 

52.0 

52.0 .110£*03 

.827E+09 

.998EKW 

50 

1 

CHUT* 

.40I48E-02 

.31613£*01 




T 

1(1) 

Y2B 

DHTH 

R2 

N 

CONTROL LOOP 

57.0 

57.0 . I03E+C3 

- 122E*10 

•998E*00 

55 

1 

DMT* 

.40442E-02 

. 31I18E+01 




T 

1(1) 

Y2B 

DHTH 

R2 

N 

CONTROL LOOP 

02.0 

62.0 .962E+02 

•173E+10 

.998E*00 

60 

1 

CHAT= 

.40541E-02 

•31024E+01 




T 

u\> 

Y2B 

DHTH 

R2 

N 

CONTROL LOOP 

67.0 

67.0 .899E+02 

.240E+10 

.998E*00 

65 

1 

CHAT= 

.41209E-02 

. 29752EK»1 





INPUT FILE ? 


d = 0 . 001 H 21 
c = 2.975 
k = O. 19 I 4.9 
A = 722. 
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GENERIC P-8 DIGITAL SIIOLATION 


In Examples A and B of Section IV, the equations of motion for a 
generic F-8 aircraft with known "pilot" control laws were used to generate 
the time histories. The time histories were then used as the input to the 
NIPIP identification scheme. These equations of motion and known control 
laws were simulated on a digital computer. The equations of motion were 
the same for both of these examples; however, the "pilot" control laws 
were different in the two examples. The pilot control laws which were 
used are defined in Section IV. 

The trim (i.e., straight and level) flight condition was specified at 
an altitude of 7000 ft, with an airspeed of 465 fps, giving a trim 
elevator control position of -6.1 deg and a trim pitch attitude of 
5.1 deg. The airspeed was held constant at 465 fps throughout the time 
history by eliminating the axial acceleration in the equations of 
motion. The stability derivatives of this generic F-8 aircraft are listed 
In Table A-l. The hare airframe short period natural frequency and 
damping ft * this flight condition are 3.48 rad/sec and 0.6, respectively. 

In this digital simulation the aircraft was assumed to be flying 
through atmospheric turbulence which was modeled by passing a uniform 
distribution through a first-order filter. That is, 


w 

JL 

n 


o /l"7u 

w WO 

s + 2D /L 
o w 


where Wg is the vertical turbulence and n Is a random number with a 
uniform distribution. The scale length, L^, was set at a constant 
1 000 ft, an«j the rras level, a w , was 3.0 fps. 
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TABLE A-l 

LONGITUDINAL STABILITY AND CONTROL 
DERIVATIVES FOR THE F-8 AIRCRAFT 


*u - 0 2 „ 

*» - 0 z » 

* q ■ 0 2 , 

*« e ■ 0 


-0.0578 (rad/sec) 

-0.996 (rad/sec) ^ 
0 ^ 
-63.325 (fps^/rad) 


0.002 (l/sec-ft) 
-0.0204 (l/sec-ft) 
-3.187 (1/sec) 

-7.451 ( rad / sec ^/ rad) 


TRIM CONDITION: 

U Q * 465 fps 

y Q * 0.0 deg 

0 O * 5.1 deg 

» -6.1 deg 

e o 


A-2 
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OPERATING PROCEDURES 


This appendix presents a description of the disk and file operations 
necessary to run the NAS software. This description will be presented for 
two computer systems. The first is concerned with operations on an 
Apple II-Plus micro-computer and the second is for operation on the Dryden 
Cyber computer. 


A. OPERATION ON AN APPLE II-PLUS 

KICRO-ODirUTER 

A prerequisite for using the NAS software on an Apple micro-computer 
is that the Apple must have 64K of random access memory (RAM) and must 
have the Microsoft Softcard CP/M operating system and Microsoft 
FORTRAN-80. This last requirement is solely due to the fact that this is 
the operating system and FORTRAN language used to write this particular 
version of the NAS software. However, Microsoft FORTRAN-80 is fairly 
standard FORTRAN, hence with minor modifications it may be possible to use 
this software on any computer with a standard FORTRAN language capability. 

For additional information about the use of the CP/M text editor, 
FORTRAN-80 compiler and linker, or other operating system features, refer 
to Refs. B-l through B-3. 


1. The NAS Program Diskette 

The NAS programs reside on one diskette. The files on the NAS 
diskette are: 
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a. FORTRAN sources files (.FOR files) 
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NSAP.FOR — Main line program 
NIPIP.FOR — Subroutines 
FUNC . FOR — Functions 


MODE. FOR 
PREPR . FOR 


Default versions of the user definable subroutines 


b. Relocatable Object Files (.REL Files) 
NSAP.REL 
NIPIP.REL 
FUNC. REL 
MODE. REL 
PREPR. REL 


c. Executable Object Files (.COM Files) 
NAS.COM 


d. Submit Files (.SUB Files) 
FL80.SUB 


The FORTRAN source files (.FOR) consist of programs as typed in by the 
user. (Listings of these programs, in the form of microfiche, can be 
found in the envelope inside the back cover of this report.) The 
relocatable object files (.REL) are machine readable binary files produced 
by the FORTRAN-80 compiler from the FORTRAN source programs. The 
executable object file NAS.COM is the absolute version of the NAS program 
that Is properly loaded and ready to be executed. This file is the only 
file which is necessary for program operation; however, the .REL and .FOR 
files are required to arrive at this executable version. 
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FL80.SUB is a program which defines and executes macro instructions. 
FL80.SUB compiles and links the NAS software on a multi-drive system, 

i.e., performs the steps necessary to convert the source files into the 
executable version of the software. 


2. Hiking a Backup Diskette 

Always work with a backup copy of the NAS system master and store the 
original for safe keeping. Thus there will always be a default version of 
the software available. The procedure for making a backup is outlined in 
Ref. B-3. The following description relies on the fact that the backup 
copy is being used, as some files need to be erased when compiling and 
linking edited versions. 


3. Sample Session for a Single-Drive System 

a. Organization of Diskettes , uue to lack of disk space on a single 
diskette, the NAS programs and the FOFTRAN-80 software will not fit on a 
single diskette. Thus it is necessary to use two diskettes. One contains 
a copy of the NAS programs listed above as well as the CP/M text editor. 
The other contains the FORTRAN-80 software. For convenience a third 
diskette contains the input data file. 


NAS Disk 


FORTRAN-80 Disk 


.FOR Files 
.REL Files 
.COM Files 
.SUB Files 
ED.COM 


F80.COM 

L80.COM 

FORLIB.REL 
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b. Editing the Source Programs * The NAS program contains two user 
def*~able subroutines, MODE and PREPR, and at some point it may be 
necessary to change these programs from their default versions. 

It should be noted; although it is possible to edit any of the source 
programs, it is not recommended or required to edit any source programs 
other than MODE or PREPR for program operation. The source files are 
edited using the CP/M editor, ED*, by typing 

ED FN.FOR 

where FN is the file name of the source file to be edited, i.e., MODE or 
PREPR. 

After the source file, FN, has been edited, exit the editor by typing 


E 


c. Compiling and Linking Edited Source Programs . The next step is to 
compile and link the edited program. This procedure is explained and 
outlined in the FORTRAN-80 user's guide, Ref. B-l. The commands and the 
procedure are summarized below: 

Place the FORTRAN-80 disk into the disk drive and type 

F80 

When the asterick prompt appears, place the copied NAS disk into the drive 
and type 


«* FN 


*Refer to Ref. B-3 for text editor commands. 
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where FN refers to the name of the source file to be compiled, i.e., PREPR 
or MODE. This creates the .REL object code with the same name. Then type 
control-c and enter* 


ERA NAS.COM 

Place the FORTRAN-80 disk into the drive, type control-c, and enter 

L80 

When the asterisk appears, Insert the NAS disk and enter 
NSAP .NIPIP , FUNC, MODE , PREPR , FILENAME/N 

where filename specifies the name to be given to the executable object 
file. Note if future editing is to be done on this diskette and a 
filename other than the one previously used is specified, it is necessary 
to erase the earlier version of the executable object code as was done 
when NAS.COM was erased. This is due to disk space limitations. Thus 
substitute 


ERA FILENAME.COM 

for ERA NAS.COM in this procedure and specify the new file name in the 
line after the L80 command. 

When finished, the linker will return an asterisk. Do not worry about 
the list of "undefined globals." Put the FORTRAN-80 disk back in the 
drive and enter 


*First make sure this is a back-up copy of the master diskette, 
because this command will erase the default executable NAS program. If it 
is not a backup, first create a backup copy of the master diskette. See 
Ref. B-3 for details on how to do this. 
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FORLIB/S 


This command searches the FORTRAN library and links the program with the 
required subroutines. When the asterisk prompt appears, insert the copied 
NAS disk into the drive and type 


/ E 

This will save the machine-executable file Filename.COM on the disk. 

d. Run the Program . After the program has been compiled and linked, 
it is executed by typing 


FILENAME 


When the Apple responds 


INPUT FILE? 


and the red light goes out on the disk drive, remove the copied NAS disk 
and place the diskette containing the input data file into the disk 
drive. Then enter the filename of the input data file and the program 
execution continues as shown in the examples of Section IV. 

This ends the sample session for a single-drive system. The next 
section presents a similar discussion for multi-drive systems. 


4. Sample Session on a Hilti-Drive System 

a. Organization of Diskettes . Due to disk space limitations, two 
diskettes are required for the FORTRAN-80 programs and the NAS programs. 
One contains a copy of the NAS programs. The other contains the 

FORTRAN-80 software, the CP/M text editor, and the file SUBMIT.COM. A 
third diskette contains che input data file. 
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NAS Disk FORTRAN-80 Disk 

(Drive A:) (Drive B:) 

.FOR Files F80.COM 

.REL Files L80.C0M 

.COM Files FORLIB.REC 

.SUB Files ED.COM 

SUBMIT.COM 

b. Editing the Source Programs . The NAS program con^xus two user 
definable subroutines, MODE and PREPR, and at some point it may be 
necessary to change these programs from their default versions. It should 
be noted, although it is possible to edit any of the source programs, it 
is not recommended or required to edit any source program other than MODE 
or PREPR for program operation. The source programs are edited using the 
CP/M text editor, ED, by typing 

B:ED FN.FOR 

where FN is the file name of the source file, i.e., MODE or PREPR. After 
the source file FN has been edited, to exit the editor type 

E 

c. Compiling and Linking Edited Source Programs . The next step is to 
compile and link the edited program. This procedure uses the macro 
command file FL80.SUB. The macro commands are F”mraarized below: 
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Place the copied NAS disk into Drive A: , the FORTRAN-80 disk Into 
Drive B:, type control-c, and type* 

ERA NAS.COM 

Bs SUBMIT FL80 FILENAME 

where FILENAME is the name to be given to the machine-executable object 
file. Note, if changing the NAS programs from their default version, it 
is recommended that some name other than NAS be used for FILENAME in the 
above command. Thus NAS. COM will always be the default version of the 
software. 

The Apple will respond 

ERA *.BAK — Erases all .BAR files created by the CP/M text 
editor. 

B:F80 ■ MODE — Compiles the MODE. FOR program. 

F:F80 - PREPR — Compiles the PREPR.FOR program 
B: — go to Drive B: 

L80 A:NSAP,A:NIPIP,A:FUNC,A: MODE.AsPREPR, A:FILENAME/N/E 

Creates the machine-executable object file. This file — called 
FILENAME.COM — is saved on Drive A: 

A: — Return to Drive A: 


* First, is this a back-up copy of the master diskette, as this command 
will erase the executable NAS program. If not, first create a back-up 
copy of the master diskette, see Ref. B-3 for details. 
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Note that if future editing is to be done on t> is diskette and a file- 
name other than the one previously used Is specified in the submit 
command , it is necessary to erase that version of the executable object 
code as was done when NAS.COM was erased. This is due to disk space lim- 
itations. Thus, substitute 


ERA FILENAME .COM 


for ERA NAS.COM in this procedure and specify the new filename in the 
submit command. 

d. Run the Program . Turn on printer if hard copy is desired. After 
the program has been compiled and linked, it is executed by placing the 
diskette with the new executable object code into Drive B: and the input 
data disk into Drive A:, typing control-c, and entering 

B: FILENAME 


When the program responds 


INPUT FILE? 


Enter the input data file name and program execution continues as shown in 
the examples of Section IV. 


B. OPERATION ON THE DRYDEN CYBER COMPUTER 

The NAS software resides on the Cyber in two forms: one is an inter- 
active version and the other a batch version. The only difference between 
these two versions Is that, In the batch version, all NAS variables must 
be defined prior to execution using subroutine MDDE. The control flag MDD 
is inoperable in the batch version, and the program always calls subrou- 
tine MDDE; whereas, In the interactive version, the user has the option of 
using the namelist or * jbroutir j MODE for defining the NAS variables. 
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Both of the Cyber versions (interactive and batch) are written in 
standard FORTRAN hence should run on any computer with FORTRAN capability 
with minor modifications. FORTRAN-V is the Cyber standard FORTRAN (see 
Ref. B-4); however, the NAS software does not rely on any special capabil- 
ities of FORTRAN-V, hence the NAS software can be run using FORTRAN IV. 

1. The HAS Program Files 

The two NAS versions reside in two separate files. The source files 
on the Cyber are: 


NASIAF — Interactive FORTRAN version 
NASSUB — Batch version 


These files are the FORTRAN source files, and they contain the mainline 
program, all function subprograms, subroutines, and the default user- 
defined subroutines, MODE and PREPR. Thus the NAS program is self- 
contained in the above files for both applications. 


2. Miking a Backup Source File 

Always work with a backup copy of the NAS system master and store the 
original for safe keeping. Thus there will always be a default version of 
the software available. The procedure for making a backup is outlined in 
Ref. B-5. The following description relies on the fact that the backup 
copy is being used, as some sections may need to be edited when adapting 
the program to the problem at hand. 
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3. Staple Session Using the Interactive Version 


a. Calling the Interactive Version . The following commands assume 
that you have a copy of the NASIAF version on your disk. If not, obtain a 
suitable copy. After logging onto the Cyber (see Ref. B-5). The Interac- 
tive FORTRAN (IAF) version is run using the FORTRAN compiler. 

Enter the FORTRAN compiler by typing 

FTNTS. 

IAF responds with the following: 

OLD, NEW, OR LIB FILE: 


Enter OLD, NASIAF 
IAF responds : 


READY. 


The ready response Indicates that IAF is ready to accept another 
command . 

b. Editing the Source Programs . The NASIAF program contains two 
user-definable subroutines, MODE and PREPR, and at some point is may be 
necessary to change these programs from their default versions. 

It should be noted, although it Is possible to edit any of the source 
programs, it Is not recommended or required *•« edit any source programs 
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other than MODE or PREPR for program operation. The source file is edited 
using XEDIT* by typing 


XEDIT, NASIAF 

where NASIAF is the file name of the Interactive source file to be 
edited. The user-defined subroutines are at the end (bottom) o r the 
source file. They can be found using the locate command. 

After the source file has been edited, exit the editor by typing 

Q, FILENAW, RL 

where FILENA1C Is the filename to be given to the edited source file. RL 
will replace the file and make the file a local file. 

c. Attaching the Data Input File . The next step is to attach the 

data input file. The attach procedure is explained and outlined In the 
NOS Reference fenual, Ref. B-5. The commands and the procedure are 
summarized below: 

The default Input file Is DATA; thus input data files are required 
to be attached under the name DATA using 

ATTACH, DATA - INFILE 

Where INFILE is the file name of the input data file to be 
attached . 

d. Run the Program . After the input data file has been attached, the 
program Is executed by typing 


RUN 


*Refer to Ref. 


B-6 for text editor commands. 
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When the Cyber responds 


INPUT FILE? 

Type the name of the input file (i.e., INFILE) This is used only for 
future reference when reviewing the program output; the program will use 
the file attached as DATA. Program execution continues as shown in the 
examples of Section IV. 

This ends the sample session for the interactive version. next 

section presents a similar discussion for the batch version. 


4. Sample Session Using the Batch Version 

a. Calling the Batch Version . The following commands assume that you 
have a copy of the NASSUB version on your disk. If not, obtain a suitable 
copy. The batch version is run using an interactive active terminal and 
submitting the program as a batch job. Thus all NAS variables must be 
defined prior to program execution using subroutine MODE. After logging 
onto the Cyber, the batch version is obtained by typing 

FTNTS 


IAF responds 


OLD, NEW, OR LIB FILE; 


Enter OLD, NASSUB 
IAF responds 


READY. 
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The ready response indicates that IAF is ready to accept another 
command • 

b. Editing the Source Programs . Since the batch version requires the 
NAS variables to be stored in subroutine MODE, it is always necessary to 
edit this subroutine to define the appropriate NAS variables. It may also 
\ necessary to change subroutine PREPR from its default version. It 
should be noted, although it is possible to edit any of the source pro- 
grams, it is not recommended or required to edit any source programs other 
than MODE or PREPR for porgram operation. The source file is edited using 
XEDIT* by typing 


XEDIT,NASSUB 

where NASSUB is the file name of the interactive source file to be 
edited. The user-defined subroutines are at the end (bottom) of the 
source file. They can be found using the locate command. 

After the source file has been edited, exit the editor by typing 

Q, FILENAME, RL 

where FILENAME is the file name to be given to the edited source file. RL 
will replace the file and make the file a local file. 


* 


Refer to Ref. 


B-6 for text editor commands. 
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c. Batch Job Control Qards . The batch version, NASSUB, contains job 
control cards. These cards are summarized below: 


/JOB 



JOB NAME. 

Name to be given to 

Job 

USER, USERNAME, PASSWORD. 



CHARGE, N01 , N01 , UTILITY. 

Charge numbers 


SETTL, 2000. 

Sets time limits 


FTN5(ET«F,PL,L«0). 

Use the FORTRAN V compiler 

ATTACH, DATA - INFILE. 

Attaches input data 

file 


INFILE under name 

DATA 

LGO. 

Load and execute 


/EOR 



NASSUB PROGRAM 



/EOF 




These cards must be adapted to a valid user when running the program. 

d. Running the Batch Program . After the proper job control cards 
have been specified, the itch version is run by typing 

SUBMIT, FILENAME, B 

where FILENAME is the name given to the program when exiting the editor, 
and B specifies the batch mode. The default specification for the program 
output device is the line printer. If any other device is to be speci- 
fied, the control flag IOUT needs to be changed as well as specifying the 
suitable device in the program card in the mainline program 

This ends the sample session for the batch version. 
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APPENDIX C 


SUMMIT OF HAS PROGRAM VARIABLES 


Table C-l contains summary descriptions of all of the NAS param- 
eters. More complete descriptions of the parameters can be found in 
Section III of the main body of this report. 

The list contained in Table C-l will be very useful to the experienced 
NAS user for a number of reasons. First, the table can be used as a check 
list for specifying the NAS parameter values. Second, the "Variable No." 
column in the list can be used when inputting the values using the name- 
list. Finally, the list significantly speeds up the time required to use 
the NAS program. 
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SUMMARY OF NAS PROGRAM VARIABLES 


Variable 

Name 

Variable 

No. 

Default 

Value 

Frequently 

Changed? 

Brief Description 

Page 

No. 

IN 

18 

10 

No 

Logical unit for input device 

26 

IOUT 

19 

6 

No 

Logical unit for output device 

26 

MOD 

20 

0 

No 

"Director” or program control 
for inputs of NAS variables 

26 

NPS 

21 

1 

No 

Number of control loops 

26 

INA 

22 

1 

Yes 

Verify F-vector 

27 

ISA 

23 

1 

Yes 

Compute frequency response 

27 

ISKIP 

24 

1 

No 

Debug option 

28 

INFW 

25 

1 

Yes 

Sliding time window ON (0), 
OFF (1) 

28 

TST 

1 

0.0 

Yes 

Start time 

29 

TEND 

2 

100.0 

Yes 

Stop time 

29 

TDT 

3 

0.1 

Yes 

Frame time 

30 

TOUT 

4 

5.0 

Yes 

Output time 

30 

TC 

5 

5.0 

Yes 

Compute time 

30 

TW 

6 

— 

Yes 

Time window 

30 

NC 

26 



Yes 

Number of DOB’ 

32 

IBIAS 

27 

— 

Yes 

Bias DOF, ON (1), OFF (0) 

32 

NVAR 

28 

— 

Yes 

Number of independent variables 

32 

INDEX j 

index 2 

29 

— 

Yes 


32 

30 

— 

Yes 

Indirect address of var- 

32 

INDEX 3 

31 

— 

Yes 

iables used in the estimation 

32 

index 4 

32 

— 

Yes 

equations. Length is NVAR 

32 

index 5 

33 

“ 

Yes 


32 

NVECj 

nvec 2 

NVEC 3 

nvec 4 

34 

— 

Yes 


33 

35 

— 

Yes 


33 

36 

— 

Yes 


33 

37 

— 

Yes 

Specify format of the 

33 

nvec 5 

38 

— 

Yes 

F-vector. Length is NC 

33 

nvec 6 

39 

— 

Yes 

if IBIAS - 0, NC-1 if 

33 

nvec 7 

40 

— 

Yes 

IBIAS - 1. 

33 

NVECg 

41 

— 

Yes 


33 

NVEC 9 

‘JVECjq 

42 

— 

Yes 


33 

43 

” 

Yes 


33 

BV! 

13 



Yes 


34 

bv 2 

14 

— 

Yes 

Reference values of the var- 

34 

bv 3 

15 

— 

Yes 

iables in the estimation 

34 

bv 4 

16 

— 

Yes 

equation. Length is NVAR 

34 

bv 5 

17 

— 

Yes 


34 
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Variable 

Variable 

Default 

Frequently 


Page 

Name 

No. 

Value 

Oianged? 

Brief Description 

No. 

WMIN 

7 

0.1 

No 


35 

WMAX 

8 

10.0 

No 


35 

WINC 

9 

2.0 

No 


35 

HFG 

10 

1.0 

No 


35 

NB 

44 

— 

No 


35 

NA 

45 

— 

No 


35 

KTAU 

46 

— 

No 


35 

NBI, 

nbi 2 

47 

— 

No 

Parameters used to compute 

35 

48 

— 

No 

frequency response 

35 

NBI 3 

49 

-- 

No 


35 

nbi 4 

50 

— 

No 


35 

NBIc 

51 

— 

No 


35 

NAIj 

nai 2 

52 

— 

No 


35 

53 

— 

No 


35 

NAI, 

54 

— 

No 


35 

NAI 4 

55 

— 

No 


35 

nai 5 

56 

— 

No 


35 

IFOPT 

57 

— 

No 


37 

TOL 

11 

1000.0 

No 

Warning tolerances 

37 

RLD 

2 

0.0001 

No 

Warning tolerances 

38 
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